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ELF/VLF/LF  PROPAGATION  AND  SYSTEM  DESIGN 


INTRODUCTION 

The  utilization  of  the  longwave  channel  (5  Hz  to  300  kHz)  for  digital  radio  transmission  has 
greatly  increased  during  the  past  20  years.  Several  large" VLF  ground  based  transmitting  antennas  have 
been  rebuilt,  and  new  facilities  have  been  constructed.  The  ELF  band  has  been  opened  up  for  use. 
Squadrons  of  longwave  broadcasting  aircraft  have  been  procured  and  operated.  A  widespread  use  of 
new  modulation  techniques  and  receivers  has  been  verified.  The  importance  of  systems  using  the 
longwave  band  is  not  widely  appreciated,  but  the  statement  by  Donald  C.  Latham,  Assistant  Secretary 
of  Defense  for  Command,  Control,  Communication  and  Intelligence  in  the  Appendix  (1],  is  abundantly 
clear.  It  is  impossible  to  present  a  comprehensive  picture  of  the  long  wavelength  situation  today 
without  producing  a  book  as  large  as  A.  D.  Watt’s  monumental  treatise  (21-  However,  because  many 
students  and  readers  are  unfamiliar  with  the  longwave  area,  we  present  a  brief  introductory  discussion 
of  the  total  longwave  communications  system  including  transmitting  antennas  and  radio  receivers 
before  discussing  in  detail  the  impact  of  propagation  on  system  design  and  performance. 

ANTENNAS 

The  longwave  channel  antennas  are  an  important  component  of  the  overall  propagation  medium 
because  they  provide  a  basic  restriction  on  types  of  signals  that  can  be  launched.  The  longwave  channel 
is  characterized  by  the  fact  that  efficient  radiation  from  conventional  size  man-made  structures  is  diffi¬ 
cult  to  achieve.  Great  ingenuity  has  been  employed  to  produce  adequately  large  antenna  configurations 
and  to  obtain  the  greatest  possible  efficiency  from  antenna  systems  that  are  necessarily  small  in  com¬ 
parison  to  a  wavelength.  Figure  1  shows  the  types  of  longwave  antenna  commonly  used. 
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Many  of  the  fixed/ground  based  vertical  monopole  antennas  are  of  the  “small  antenna"  variety 
characterized  by  a  very  small  radiation  resistance  (less  than  1  O).  Consequently,  these  antennas  are 
quite  inefficient  unless  special  precautions  are  taken  to  reduce  utterly  all  unnecessary  losses.  One  com¬ 
mon  technique  is  the  construction  of  a  large  ground  plane  to  suppress  losses  in  the  earth  near  the 
antenna.  A  second  feature  is  the  deployment  of  an  elaborate  “top  hat”  consisting  of  cables  connected 
together  and  suspended  from  towers  to  provide  additional  antenna  capacity  to  increase  wirease  the 
antenna’s  ability  to  carry  current  for  a  given  voltage.  In  this  way  the  losses  are  reduced  to  a  minimum 
and  a  higher  efficiency  is  achieved.  However,  a  high  “Q"  antenna  system  with  a  consequent  low 
bandwidth  results  from  these  procedures.  The  design  considerations  for  such  antennas  are  summarized 
in  Refs.  1  through  3.  Figures  2  through  5  (3]  show  the  four  principal  ground  based  antenna  types:  Tri- 
atic,  Goliath,  Trideco,  and  valley  span.  Figures  6  through  11  illustrate  the  designs  of  U  S.  Navy  VLF 
fixed  ground-based  antennas. 

The  suspension  of  a  long  conducting  cable  from  an  aircraft  or  other  lifting  vehicle  is  an  alternative 
method  for  making  an  efficient  radiator  at  long  wavelengths.  Thus  a  half  wave  dipole  antenna  is 
achieved  with  a  large  radiation  resistance  and  high  radiation  efficiency  without  the  high  “Q”  and  lim¬ 
ited  bandwidth  of  the  ground-based  tower  antennas.  This  technology  is  the  basis  of  the  TACAMO  and 
airborne  command  post  longwave  communication  capability  as  presented  in  Refs.  5,  6.  and  7.  The  air¬ 
borne  facilities  are  illustrated  in  Fig,  12,  [5] . 

An  horizontal  electric  antenna  can  be  used  at  long  waves.  This  antenna  is  utilized  for  the  ELF 
communication  system  (8,9] .  In  this  system,  a  powerful  current  is  driven  along  horizontal  wires  that 
like  no'i,er  lines  are  stretched  over  the  earth’s  surface  and  grounded  at  each  end.  When  the  conduc¬ 
tivity  of  the  underlying  earth  is  low,  it  is  as  if  the  antenna  is  suspended  at  a  modest  height  above  the 
“elect! ka!”  earth.  The  effective  image  currents  flow  from  one  terminal  to  the  other  at  a  depth  propor¬ 
tional  to  the  electromagnetic  skin  depth  of  the  underlying  material  This  mechanism  permits  the 
horizontal  antennas  to  operate  with  modest  efficiency. 

In  the  near  future,  it  will  be  feasible  to  deploy  a  large  longwave  antenna  from  an  earth  orbiting 
satellite.  Italy  and  the  United  States  (NASA)  in  a  joint  program  are  developing  a  60  km  conducting 
tether  system  for  use  in  space.  One  of  the  applications  of  this  tether  is  the  generation  of  long  radio 
waves  (101. 
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Fig  2  —  Trialic  type  antenna 
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I'ig  12  —  Airborne  submarine  communications  (Figure  3  on  p  28.  Vol.  39,  No.  T  entitled  "Strategic 
Command  and  Control  and  National  Security,"  by  Dr.  Bruec  G  Blair  from  Signal,  published  by  The 
Armed  Forces  Communications  and  Electronics  Association  Copyright  (c)  1985,  used  b>  permission.) 


PROPAGATION 

Summary 

The  first  transatlantic  radio  propagation  experiments  of  Marcon.  were  conducted  in  the  long- 
wavelength  band  [2]  In  comparison  to  the  short-wave  band,  the  long-wavelength  hand  is  highly  reli¬ 
able.  Long-distance  propagation  of  a  long  radio  wave  relies  on  a  reflection  from  the  lower  ionosphere 
to  bring  the  radio  wave  back  to  earth.  On  the  other  hand,  short-waves  are  gradually  bent  back  to  earth 
because  of  a  refraction  process  that  depends  on  the  detailed  horizontal  and  vertical  gradients  of  an 
upper-ionosphere  F  layer.  A  VLF  wave  spends  most  of  its  journey  in  the  free  space  between  the 
earth’s  surface  and  the  lower  ionosphere.  A  short-wave  propagates  mostly  through  the  ionosphere 
between  70  and  300  km  height,  between  the  D-laycr  and  the  F-laycr  This  is  probably  the  basic  reason 
that  long-wavelength  waves  show  greater  stability  of  reception.  The  disadvantages  of  long  waves 
(which  were  deadly  from  an  ordinary  economic  point  of  view)  are: 

(a)  Large  antenna  structure 

(b1  Narrow  bandwidlhs 

(c)  Limited  total  spectrum 

(d)  High  powers  required  to  overcome  noise. 

The  advantages  of  long  waves  are; 

(a)  Great  propagation  range  (even  to  the  antipodes) 

( b >  High  phase  stability 

(c)  Significant  ability  to  penetrate  earth  and  seawater 
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The  long  propagation  range  and  phase  stability  for  long  wavelength  waves  make  their  use  for  navi¬ 
gation  feasible  (Omega,  Decca,  and  Loran).  Their  penetration  into  the  earth  makes  the  waves  useful  to 
prospectors  looking  for  ore-bodies  or  other  deposits  beneath  the  earth. 


The  Ionosphere  and  Reflection  Coefficients 


The  layers  ascribed  to  the  ionosphere  and  atmosphere  are  illustrated  in  Fig.  13  [11].  The  ionos¬ 
pheric  reflection  coefficients  are  defined  as  the  ratio  of  downcoming  electric  field  to  the  upgoing  field. 
Because  there  are  two  orthogonal  upgoing  and  downcoming  field  components,  four  complex  reflection 
coefficients  characterize  the  process  at  any  given  angle  of  incidence.  The  matrix  form  is  as  follows: 


ii^ii  i  ^ ii 
II^L  L^l 


(1) 


where 

E"  is  the  upgoing  electric  field  in  the  plane  of  incidence; 

£“  is  the  upgoing  electric  field  perpendicular  to  the  plane  of  incidence 
(i.e.,  the  horizontal  electric  field); 

Ef,  is  the  downcoming  electric  field  in  the  plane  of  incidence; 

Ef  is  the  downcoming  electric  field  perpendicular  to  the  plane  of  incidence. 
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!ig.  13  —  Atmospheric  layers  and  nomenclature 
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The  reflection  coefficients  are  functions  of  the  electron  and  ion  density  profiles  of  the  ionosphere. 
The  collision  frequency  profiles,  geomagnetic  field  and  angle  of  incidence  must  also  be  specified. 
Several  full-wave  techniques  can  be  used  to  calculate  the  reflection  coefficient  matrix.  It  is  beyond  the 
range  of  this  lecture  to  describe  these  methods  in  depth.  Figures  14  through  16  show  sample  calcula¬ 
tions  of  reflection  coefficients  versus  incidence  angle  at  3  kHz  for  a  geomagnetic  dip  angle  of  60°.  (Ref. 
12)  Figure  17  shows  two  reflection  coefficients  as  a  function  of  frequency  for  fixed  incidence  angle  of 
approximately  30°.  (Ref.  13)  The  ionospheric  profile  used  for  each  calculation  is  shown  in  Fig.  18. 
Notice  that  the  calculated  reflection  coefficients  above  about  30  kHz  are  strongly  varying  with  fre¬ 
quency.  They  seem  to  reflect  the  detailed  structure  of  the  ionosphere.  The  smooth  approximation  by 
Allcock  and  Belrose  (Ref.  14)  seems  adequate  for  a  first  approximation;  but  it  lacks  the  oscillations 
noted  in  the  fuliwave  results.  During  nighttime  long  radio  waves  are  transmitted  through  the  lower 
ionosphere  with  small  loss.  Space-borne  receivers  can  readily  hear  ground-based  transmitters.  Likewise 
receivers  on  the  ground  can  sense  many  curious  electrical  noises  that  are  generated  and  stimulated  in 
the  ionosphere  and  magnetosphere.  For  these  ground-to-space  and  space-to-ground  paths  the  lower 
ionosphere  transmission  coefficient  is  significant  in  forming  a  signal  power  budget.  This  transmission 
coefficient  can  be  calculated  using  Pitteway’s  full  wave  method.  (Ref.  15).  Sample  transmission  coeffi¬ 
cient  results  for  a  3  kHz  case  are  shown  in  Figs.  19  and  20. 

Field  Strength  Calculations  using  Waveguide  Modes 

For  long  range  propagation  calculations  it  is  customary  [16]  to  treat  the  space  between  the  earth 
and  the  ionosphere  as  a  waveguide  and  to  calculate  the  field  strengths  from  a  transmitter  as  a  summa¬ 
tion  of  waveguide  modes.  A  popular  method  for  achieving  this  is  by  using  the  WAVEGUIDE  com¬ 
puter  program. 
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~ig.  IS  —  3  kHz  nighttime  ionospheric  reflection 
coefficients,  south-to-north  propagation 


Fig.  16  —  3  kHz  nighttime 

ionospheric  reflection  coefficients, 
easl-lo-wesl  piopagation 
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Fig  17  —  Reflection  coefficients  vs  frequency  The  curve  is  a  plot  if  an 
empirical  reflection  coefficient  approximation  given  by  Belrose  (14)  using  an 
equivalent  frequency  concept  of  Allcock-  The  n/?1and1/?1  reflection 
coefficients  follow  the  convention  of  Ref.  IS. 
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Fig.  18  —  Electron  density  profile 
for  nighttime  calculations 
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Description  of  WA  V EG  VIDE 

The  WAVEGUIDE  computer  program  originally  was  developed  at  the  Naval  Electronics  Labora¬ 
tory  Center  to  predict  the  vertical  electric  field  received  at  a  point  on  the  earth's  surface  and  produced 
by  a  vertical  electric  transmitting  antenna  at  another  point  on  the  earth's  surface.  Later  it  was  extended 
to  give  the  crosspolarized  (horizontal)  components  of  the  field  at  any  height  in  the  earth-ionosphere 
waveguide  and  even  in  the  ionosphere  (17-19],  The  basic  formulation  of  the  equations  for  the  fields  is 
given  in  Ref  20.  The  net  resultant  vertical  electric  field  at  the  earth's  surface  is  represented  as  a  sum¬ 
mation  of  waveguide-mode  fields. 


E. (0)  = 


sin  (d/a) 


Lisin1'2  8„  A „  exp  (in/ 4  +  /X0r/sin  0,)], 


where 


8„  is  the  eigenangle  of  the  r?th  waveguide  mode. 

Ids  is  the  dipole  moment  (ampere-meter)  of  the  vertical  antenna  located  on  the  earth's  sur¬ 

face, 

A,  is  the  excitation  factor  of  the  nth  waveguide  mode, 

a  is  the  radius  of  the  earth, 

d  is  the  great-circle  distance  between  the  transmitter  and  the  receiver, 

h  is  the  reference  height  of  the  ionosphere  (used  in  the  definition  of  \„). 

kn  (=  co/c  )  is  the  free-space  propagation  constant  of  the  wave, 

X  is  the  wavelength  of  the  wave, 

and 

p.u  and  are  the  magnetic  permeability  and  dielectric  constant  of  free  space. 

Each  waveguide  mode  is  characterized  by  its  own  value  of  attenuation  rate  «„  and  phase  velocity 
v;.  which  arc  related  to  the  eigenangle  of  the  nth  waveguide  mode  h„  according  to 

a„  =  0.02895  w  !m  (sin  8,,)  (3) 

and 

v"  =  - - -  (4) 

p  Re(sin  8„) 

where  «„  is  measured  in  decibels  per  1000  km  (dB/Mm)  and  c  is  the  speed  of  light 


The  values  of  the  waveguide-mode  eigenangles  depend  on  the  reflection  coefficients  of  the 
ground  and  the  ionosphere.  Because  of  its  anisotropy,  the  ionosphere  has  four  reflection  coefficients 
,-R, ,(&„),  /(.  (8,.),  ,  /?,  (f)„),  and  R,,(8r)  for  a  given  angle  of  incidence  9„.  The  ground  is  assumed  to 
be  isotropic  and  to  have  reflection  coefficients  R  (9n)  and  ;/?,(0„).  that  are  calculable  from  the 
eigenangle,  ground  conductivity,  and  dielectric  constant.  The  WAVEGUIDE  program  uses  a  procedure 
for  calculating  the  ground  and  ionospheric  reflection  coefficients  referenced  to  any  height  z  in  the 
waveguide.  References  26  27,  and  28  discuss  the  reflection  coefficients  furthet.  The  eigenangle  8n  for 


a  given  waveguide  mode  is  obtained  by  satisfying  the  mode  equation 

If,  0  1 

1  oil 

g 

°~  If..*,  <»„> 

0 

■  |o  1  r 

(5) 

s- 

This  matrix  equation  is  simply  rewritten 

r(9„)  =  (,«,  R-.  -  I)  (;  Rl  „  -1)  -  t  R  ,  '  =  o.  (6) 

where  the  dependence  of  the  reflection  coefficient  on  8n  has  been  suppressed  to  simplify  the  notation. 
The  new  variable  E(8„)  is  also  defined,  tach  value  of  8n  which  satisfies  Eq.  (6)  is  the  eigenangle  of 
the  «lh  waveguide  mode.  The  excitation  factor  A „  for  the  nth  waveguide  mode  is  obtained  from 
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a  .  kh  .  . 

A„  -  -/  —  sin  9„ 


0  +■■*..)*  (1  -,RllRl) 

p  BF(e) 

"R"  Be 


where  the  variable  F(.0)  is  as  defined  in  Eq.  (6).  In  addition  to  the  vertical  electric  field  Et  for  a  given 
waveguide  mode  having  eigenangle  0„,  in  general,  there  exist  five  other  nonzero  field  components  £x, 
Ey,  Hx,  Hy ,  and  Ht  within  the  waveguide  for  each  mode. 

These  extra  fields  arise  because  the  ionosphere  is  anisotropic  because  of  the  earth’s  magnetic 
field.  Propagation  through  and  reflection  from  such  an  anisotropic  medium  rotates  the  plane  of  polari¬ 
zation  of  the  incident  wave  to  generate  crosspolarized  reflected  fields.  Alternatively  one  could  consider 
the  linearly  polarized  wave  transmitted  from  the  vertical  electric  dipole  antenna  to  be  a  linear  superposi¬ 
tion  of  right-handed  and  left-handed  circularly  polarized  waves,  each  reflected  with  a  different  ampli¬ 
tude  at  the  ionosphere.  The  ionospheric  reflection  causes  a  mixture  of  wave  polarizations  at  a  reception 
point  on  or  above  the  earth’s  surface.  The  fields  at  height  z  are  related  to  the  vertical  electric  fields  at 
the  earth’s  surface  £,(0)  in  a  given  waveguide  mode  by  the  following  functions  of  reflection  coeffi¬ 
cients  and  Hankel  functions: 

£,(z)/£,(0) -/„  (z),  (8) 


£,(z)/£.(0)  =  g(z)/S. 


£y(z)/Et(0) 


_ n^l  (H~i  Rj  )  A 

5(1  +  „*„)  (1  -  lR  ‘lRl)  ' 


U-M1  +i*i>/i<*> 

H,  (z)/ E,  (0)  - - - - — - = - -, 

T,(l  +„£.,)  (1  ~LRnRL) 


H„  (z)/ £  (0)  -  , 

17  S 


Hx(z)/E,(  0) 


,,/Ml  +  ;*i) 


df,  (2) 


■n  k  5(1  +  „/?„)  (I  -  LRiLRL)  dz 


where 


/,.  (z)  —  exp 


z-D  F i  M<?>  +  Ri 

a  F]  ht(qd)  +  F2  h2(qd) 


fi  (2)  - 


F}  ht(q)  +  FA  h2(q) 
F }  h\(qd)  +  Fa  h2(qd) 


g(z)  ~  —  -j-  l/„l(z)l, 

:k  dz 
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in  which 


F,  =  - 


HM  ~  i  ^ 


ak 

2 


1/3 


(A^„2  -  S2)'/2  M<?o) 


(17) 


F:  -  W,(fl0)  -  <  ~7 

N* 


f\  (N?  -  S2)'n  /i,(?0). 


1/3 


(18) 


$ 


and 


with 


3  “  “  |/l\(<?0)  ”  ' 
F 4  =  /j'i (^o>  -  » 


a/; 


1/3 


(<V2  -  S2)1/2  MW 


ak 

2 


1/3 


(N/  _  S  )  MW* 


1 

-2/3 

ak 

C2  -  -( h  -  z) 
a 


= 


2 

-2.3 

C2  -  -ih  -  D) 

ak 

a 

Q  o 


H,  (?)  -  *',<«>  +  y 


2 

-2/3 

c2  -  — 

ak 

a 

_2_ 

ak 


2  3 


h,(q)s  j  =1,2 


(19) 

(20) 

(21) 

(22) 

(23) 

(24) 


§ 

I 


Hi 

*r? 


and 


n2  =  1 - ( h  —  z ). 


«o=l - h, 

a 


K  }  I  .  (T 

a;  =  dt0  -  j  — 


(25) 


(26) 


(27) 


In  Equation  27  ir  and  e  refer  to  the  conductivity  and  dielectric  constant  of  Ihe  earth. 

Typical  Waveguide  Mode  Parameters 

Using  the  WAVEGUIDE  program,  one  can  calculate  the  propagation  parameters  of  the  various 
modes.  Figures  21  through  23  [21]  give  the  excitation  factors  and  attenuation  rates  for  the  three  lowest 
order  waveguide  modes  in  the  ELF'  band  and  lower  VLF  band.  Figures  24  through  26  show  the 
behavior  of  the  magnetic  fields  vs  height  in  the  waveguide  and  inside  the  lower  ionosphere.  The  night¬ 
time  ionosphere  as  shown  in  Fig.  18  was  used  The  direction  of  propagation  is  toward  magnetic  north 
for  a  dip  angle  of  60°.  Using  these  propagation  parameters  and  Eqs.  (7-27),  one  can  estimate  the  held 
strength  produced  by  a  transmitter.  Additional  equations  and  sample  parameters  for  this  task  are  given 
in  Ref  16 
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EXCITATION  FACTOR  UBt  ATTENUATION  RATE  (dB/Mm) 
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Figures  27  through  35  give  propagation  parameters  for  the  first  three  waveguide  modes  calculated 
for  a  typical  daytime  isotropic  ionosphere  and  used  in  the  noise  prediction  computer  program  described 
in  Ref.  22.  These  are  similar  to  but  more  comprehensive  than  the  parameters  given  in  Ref.  23, 

Field  Strength  Calculations  and  Measurements 


As  shown  in  Ref.  24,  field  strength  vs  distance  curves  are  calculated  from  these  parameters.  Fig¬ 
ures  36  through  41  (25]  illustrate  the  pattern  of  field  strength  vs  distance  for  1  kW  radiated  power 
transmitter  propagating  over  sea  water  at  24.0  kHz.  The  very  sharp  predicted  interference  hull  at  2  4 
Mm  is  noted  in  Fig.  36.  Using  the  parameters  given  in  Ref.  23,  a  statistical  prediction  of  field  strength 
vs  distance  was  made  by  a  Monte  Carlo  technique.  A  Gaussian  distribution  of  height  h  and  conduc¬ 
tivity  gradient  fi  was  assumed  to  have  an  average  h  -  70  km  and  average  0  —  .5  km-1.  The  standard 
deviation  of  h  was  1  km  and  of  $  was  .05  km"1.  This  method  predicted  the  greatest  variability  of  signal 
level  to  occur  near  the  interference  nulls  at  2.4  and  3.8  Mm.  Experimental  data  were  taken  for  a  period 
during  summer  1970  near  the  2.4  Mm  null  produced  by  the  NBA,  Balboa  transmitter  that  was  operating 
on  24.0  kFIz  (26],  The  measured  standard  deviation  was  much  smaller  than  that  predicted  from  the 
model.  It  was  also  smaller  than  that  observed  with  long-term  data-taking.  From  this  data,  it  was  con¬ 
cluded  that  during  the  short  period  of  this  experiment  the  ionosphere  was  highly  repeatable,  even 
though  it  showed  greater  variability  over  longer  durations. 

DCA  Recommended  Ionospheres 


From  an  analysis  of  field  strength  vs  distance  data  similar  to  those  shown  in  Figs.  42  and  43,  the 
US  Defense  Communications  Agency  (27)  recommended  the  ionosphere  profiles  shown  in  Table  1. 
Graphs  showing  these  ionosphere  parameters  are  shown  in  Figs.  44  and  45,  (27).  The  electron  density 
N(z)  vs  height  z  in  kilometers  for  a  particular  (i  and  h  is  given  by  Eq.  (28). 


,V(:) 


1.43  •  107-~  •  exp  (-0.15  h) 
cm1 


exp  Uj3  -  0.15)  (z  -  6)] 


(28) 


Table  1  -  DCA  [27]  Recommended  Profiles  to  Use  in  WAVEGUIDE 
OR  WAVEHOP  Propagation  Programs 


DAYTIME 

NIGHTTIME 

Summer 

Winter 

Summer 

Winter 

High  Latitude 

Low  Latitude 

fi  =  0.3,  h'  =  72 

0  -  0.5, h '  -  70 

0  -  0.3,6'  -  72 

0  _  0.3,6’  -  72 

0  =  0.5, 

6'  -  87 

The  collision  frequency  vs  height  is  given  by  Eq.  (29). 

v  (z)  =  »/0exp  (-az)  (29) 


where 


v0  =  1.82  x  l  o' •  collision/s  a  =  0.15  km  1 

From  an  analysis  of  multifrequency  midlatitude  day  and  nighttime  data  [28]  that  was  conducted  later 
during  flights  from  Hawaii  to  California,  recommendations  were  derived  which  are  given  in  Table  2  It 
is  significant  that  for  nighttime  the  laterally  homogeneous  ionospheric  profile  must  be  assumed  to 
changes  with  broadcast  frequency.  This  unphysical  situation  reflects  our  current  lack  of  understanding 
of  the  nighttime  lower  ionosphere. 
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Fig.  37  —  Histograms  of  calculated  field  strengths 
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Fig  41  —  Histograms  of  calculated  field  strengths 
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Fig.  45  —  Nighitime  Electron  Density  Profiles  and  Collision  Frequency  Profile 


Table  2  —  DCA  (28]  Recommended  Electron  Density  Profiles  for 
Use  in  Propagation  Predictions  (midlatitude) 


DAYTIME 

Winter 

NIGHTTIME  j 

Winter 

Frequency 

Profile 

Frequency 

Profile  ; 

(kHz) 

0  (km  ')  ,  H  (km) 

(kHz) 

0  (km-') 

//'  (km)  ■ 

9-60 

0  -  0  3,  //'  =  74 

below  10 

0  -  0  3. 

«■»« 7  ; 

i 

l 

10  -  !  5 

0  -■  0  4. 

H  -  87  j 

15  -  25 

0  -  0.5. 

//'-  87  : 

1 

|  25  -  30 

0  -  0.6. 

//’  “  88  ; 

30-40 

a  -  0.7. 

//'  =  88  j 

. . 

40-60 

0  -  0.8. 

(1 

<y> 

oo 

Coverage  Prediction  Programs  for  Inhomogeneous  Pains 

Up  to  this  point  we  have  discussed  prediction  methods  for  homogeneous  paths  nriiv.  However, 
the  earth's  surface,  earth's  magnetic  field,  and  the  ionosphere  typically  vary  with  distance  along  a  pro¬ 
pagation  path  References  3?  and  33  represent  early  efforts  to  apply  mode  theory  to  these  inborn  >enc- 
ous  paths  References  34  through  3'  describe  further  developments  and  elaboration  of  me  wav, guide 
mode  propagation  picture  using  mode-matching  and  the  WKB  approximation  to  propagate  the  waves 
through  spatially  varying  portions  of  ihe  path. 

Reference  38  describes  a  Very  l  ow  Frequency  Automatic  Computation  Method  (VLFACM)  com¬ 
puter  mode).  It  is  approximate  but  it  is  inexpensive  to  run  [39|.  it  makes  predictions  with  built  in 
easily-computed  values  lor  the  attenuation  and  excitation  of  the  signal.  VLFACM  only  computes  a  sin¬ 
gle  dominant  waveguide  mode  and  operates  best  from  15  to  30  kHz  The  model  parameters  were 
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selected  by  comparing  the  program  predictions  with  a  great  many  propagation  paths.  VLFACM  con¬ 
tains  expression*:  that  permit  it  to  give  a  complete  diurnal  prediction  taking  solar  zenith  angle,  geomag¬ 
netic  azimuth,  and  ground  conductivity  into  account. 

DC  A  Arctic  Propagation  Recommendation  and  Recent  Experiments 

The  DC  A  recommendations  (Ref.  27)  for  Arctic  or  high-latitude  propagation  conditions  were 
based  on  the  comparisons  shown  in  Figs.  46  through  53.  Notice  that  in  all  cases  the  field  strength  pred¬ 
ictions  are  noticeably  higher  than  the  measured  values  for  waves  propagating  across  the  Greenland  Ice 
Cap  and  that  the  conductivity  of  the  ice  cap  is  assumed  to  be  10“ 5  fl /meter  in  conformity  with  refer¬ 
ences  29  and  30.  However  at  25  kHz  reference  31  gives  the  measured  conductivity  of  an  antarctic  gla¬ 
cier  as  5.0  to  6.1  x  10~5  fi/1  m  with  a  relative  dielectric  constant  of  approximately  11.  Referring  to 
Figure  27  we  see  that  this  conductivity  change  could  raise  the  attenuation  rate  over  Greenland  at  24 
hKz  from  9dB  for  a  =  10-5  fl/m  to  between  13  and  15  dB/Mm  for  the  larger  values  of  a  given  by 
Ref.  31.  Recent  measurements  of  the  24.0  kHz  signal  from  NAA  indicate  an  attenuation  rate  over 
Greenland  in  better  agreement  with  the  higher  attenuation  rate,  suggesting  that  the  conductivity  and 
dielectric  constant  of  the  Greenland  Ice  Cap  may  need  modification. 


j _ l _ i _ i _ v _ j 

I  2  3  4  S 

DISTANCE  (Mm) 


Fig.  46  —  High-Latitude  Propagation 
Fllghi  Path:  England  to  Sondrcsirom.  Greenland  (daytime) 
Ground  conductivity:  Seawater.  4  11  /nr.  Greenland  ice  cap. 

itr!  n/m 

Transmitter:  GYN  Preston.  England  (45  kHz) 

Date:  July  1970 

Electron  Density  Profiles:  (i  =  0  3  km'1,  h'  =  72  km  with  the 
base  of  the  ionosphere  at  30  km  and. 
/3  =  0  25  km  '.  h'  =  74  km  with  the 
base  of  the  ionosphere  at  40  km 
Method  of  Computation:  WKB  approximation 


35 


AMPLITUDE.  dB  ABOVE  1  *iV/m  FOR  1  kW  RADIATED  POWER 


Fig  47  —  High-Latitude  Propagation 

Propagation  paramct(:rs  are  the  same  as  those  in  'Fig  46  except  that  only  the 
(3  =  0.3  km'1,  h'  =  72  km  profile  are  shown  Field  strengths  are  computed  for  receiver 
altitudes  (RALT)  of  0,  6.  and  9  km. 


& 


N 


KM 

*3 


S5 


s 


■  1 

f. 


ft 


> 


36 


AMPLITUDE  dB  ABOVE  1  MV/m  FOR  1 -kW  RAQIATEO  POWER 


NRL  REPORT  9028 


i  ,6  48  —  High-Latitude  Propagaiion 

I  light  Path  England  10  Sondrestrom.  Greenland  (daytime) 
Ground  conductivity  Seawater,  4  11  /mt  Greenland  ice  cap, 

I0"5  11  /m 

1  r.msmitter  MSF  Rugb>,  England  (60  kHz) 

Date:  July  1970 

I  Icctron  Density  Profiles  fi  =  0  3  km"1,  h  —  72  km  with  the 
base  of  the  ionosphere  at  30  km  and, 
ii  -  0.25  km"',  h'  =  74  km  with  the 
base  of  the  ionosphere  at  40  km 
Method  of  Computation  WKH  apprommation 
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Fig.  49  —  High-Latitude  Propagation 

Propagation  parameters  are  the  same  as  those  in  Fig.  48  except  that  only  the 
p  -  0.25  km'1,  h'  -  24  km  profile  is  shown.  Field  strengths  arc  computed  for  receiver 
altitudes  (RALT)  of  0  and  9  km. 


Fig.  50  —  High-Latitude  Propagation 
Flight  Path:  Seattle,  Washington  across  Canada 
and  over  the  Greenland  ice  cap 
ground  conductivity: 

Various  values  of  land.  sea.  and  ice  along  the  path 
Transmitter:  NLK,  Seattle.  Washington  (18  6  kHz) 
Date.  June  1958 
Electron  Density  Profiles: 
p  -  0  3  km'1,  h  -  12  km  with  the 
base  of  the  ionosphere  at  30  km  and. 

P  -  0.25  km'1,  h'  -  74  km  with  the 
base  of  the  ionosphere  at  40  km 
Method  of  Computation:  WKB  approximation 
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Fig.  51  —  High-Latitude  Propagation 
Flight  Path:  Thule,  Greenland  across  the  ice  cap 
to  London,  England 
and  Thule  to  Fairbanks,  Alaska 
Ground  conductivity: 

Various  values  of  land,  sea,  and  ice  along  the  path 
Transmitter:  NLK,  Seattle,  Washington  (18.6  kHz) 
Date:  June  1958 
Electron  Density  Profiles: 

/3  -  0.3  km'1,  h'  —  72  km  with  the 
base  of  the  ionosphere  at  30  km  and, 

/3  -  0.25  km'1,  h'  -  74  km  with  the 
base  of  the  ionosphere  at  40  km 
Method  of  Computation:  WKB  approximation 


Fig.  52  —  High-Latitude  Propagation 

The  propagation  parameters  are  the  same  as  those  of  Fig.  50  except  that  the  earth's 
conductivities  are  taken  as  only  seawater,  a  —  4.64  U/m  and  Greenland  ice,  it  ”  iO 
fl/m  The  computations  are  made  using  the  mode-conversion  procedures. 
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Fig  53  —  High-Latitude  Propagation 

The  propagation  parameters  are  the  same  as  those  of  Fig.  51  except  that  the  earth's 
conductivities  are  taken  as  only  seawater,  (  cr  “  4.64  fl/m)  und  Greenland  ice.  (  <r  ■= 
1CL5  D/m).  The  computations  are  made  using  the  mode*conversion  procedures 


Coverage  Problems 

The  question  of  the  actual  coverage  of  longwave  systems  is  currently  in  a  status  of  reevaluation 
and  flux.  Because  of  the  meager  amount  of  airborne  propagation  data  available  and  the  tendency  of  the 
ionosphere  not  to  stand  still,  we  find  that  our  coverage  prediction  techniques  are  unable  to  achieve  the 
accuracy  that  will  readily  satisfy  a  system  designer.  Significantly  large  sums  of  money  are  involved  in 
the  purchase  of  aircraft  squadrons  or  in  the  building  of  a  huge  ground  based  transmitters,  therefore 
greater  emphasis  and  resources  should  be  and  are  being  allocated  to  longwave  coverage  research. 

Wavehop  Coverage  Predictions 

The  solution  of  the  Helmholtz  equation  for  the  Hertz  vector  between  concentric  spherical  shells 
(40)  is  given  by  a  summation  of  products  of  spherical  harmonics  functions  (Eq.  30), 

jUf'  =°  f 

u  -  (2(7+1)  h™<kb)  ( kr )  Pq  (cos0)-^-.  (30) 

^  q~u  LJq 

Watson  (41)  discovered  a  technique  for  transforming  tnis  slowly  converging  series  into  an  integral  in 
the  complex  plane  and  reexpressing  the  integral  as  a  summation  over  residues. 

w  -  -/  f  •— —  /  (e)  P,[cos  it- 0)]  dv.  (31) 

J  !■  sin  err 

in  this  expression  the  sum  of  spherical  harmonic  Legendre  polynomial  Pv  (cos  (■>)  functions,  which 
utilize  integral  values  of  q,  are  reexpressed  as  an  integral  over  complex  contour  in  the  variable  v  using 
a  Legendre  function  P„  (cos  (it  -  0))  that  utilizes  a  complex  value  of  the  subscript  i>.  The  function  P, 
(cos  (rr  ~  0)]  behaves  like  a  propagating  and  exponentially  decaying  wave  in  the  <-)  direction,  i.c.. 


40 


NRL  REPORT  9028 


along  the  surface  of  the  earth.  In  the  waveguide  mode  theory  the  integral  in  Eq.  (31)  is  evaluated  by 
the  method  of  residues  at  residue  points  ut.  Each  value  of  v,  is  referred  to  as  a  waveguide  mode  root 
and  is  a  solution  of  Eq.  (32). 


D„-  0 


1  -  R™ 


R 


(v) 


( ka)hy]  ( kc ) 
( ka)hy)  (kc) 


(32) 


The  wave-hop  formulation  for  the  propagation  of  long  radio  waves  originated  at  a  very  early  lime 
after  Marconi’s  original  transatlantic  transmissions.  An  early  paper  by  Hoilingworth  [42l  is  particularly 
noteworthy  in  that  Hoilingworth  experimentally  determined  that  the  field  strength  vs  distance  of  a  long 
wavelength  station  truly  behaves  in  an  undulatory  manner,  thereby  exhibiting  the  effect  of  groundwave 
and  skywave  interference  at  ranges  governed  by  geometry  and  wavelength.  The  hop  approach  was 
adopted  by  most  workers  in  the  field  during  the  1930s  and  1940s— as  exemplified  by  Bremmcr’s  text¬ 
book  (431. 


J.  R.  Wait  (40,44]  opened  up  a  new  line  for  theoretical  advance  when  he  pointed  out  that  the 
Watson  integral  expression  for  the  fields,  (Eq.  31),  that  leads  to  the  waveguide  mode  formulation,  can 
also  be  expressed  in  terms  of  a  geometric  series  in  the  reflection  coefficients.  The  first  term  in  the 
geometric  series  corresponds  to  the  ground  wave.  The  second  term  to  a  one-hop  skywave.  The  third 
term  describes  a  two-hop  skywave,  etc.  This  concept  was  extended  by  L.  Berry,  J.R.  Johler,  and  others 
[45-49]  to  produce  a  method  for  expressing  the  fields  utilizing  the  anisotropic  ionospheric  reflection 
coefficients  in  a  realistic  way.  These  methods  were  developed  during  the  1960s  and  culminated  in  com¬ 
puter  programs  produced  by  Berry,  Herman,  and  Gonzales,  and  by  Lewis  in  the  early  1970s.  The  wave 
hop  method  has  been  shown  by  Jones  et  al.  (50]  and  Morfitl  and  Halley  (51]  to  compare  well  with  the 
waveguide  mode  methods  of  calculation. 


The  wavehop  programs  have  been  directly  applied  to  making  field  strength  prediction  maps.  The 
wavehop  method  possesses  the  great  virtue  that  after  inserting  the  parameters  of  a  propagation  path, 
the  results  arc  produced  directly  after  a  straightforward  calculation  that  requires  no  human  intervention 
to  complete.  Results  are  automatically  mapped  to  provide  signal  to  noise  coverage  charts.  Figs.  54 
through  57. 
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Fig  54  —  Coverage  churl  showing  signal  to  noise 
contours  <3  dB.  heavy  line.  0  dB  light  line)  lor  a 
naval  bioadcast  slaiion  al  Turkic.  Philippines  on 
I  19  85  kHz  lor  99%  umc  availability  for  all  hours  of 
ihe  month  ol  July.  Assumed  transmitter  power  is 
2.5  kW.  and  noise  bandwidth  is  I  kll/.. 
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Fig  57  —  Coverage  chan  for  same  conditions  as 
those  in  the  previous  Figure,  except  time 
availability  is  increased  to  99% 

Recent  Developments  m  Wave  Hop 

Recently  we  have  implemented  some  of  the  initial  ideas  of  Johler  to  make  the  wavehop  program 
amenable  to  the  treatment  of  inhomogeneous  paths  in  a  realistic  way.  To  do  this,  one  must  calculate 
ionospheric  reflection  coefficients  along  the  propagation  path  at  a  variety  of  angles  of  incidence  and 
utilize  these  values  in  evaluating  the  fields  via  interpolation  techniques.  Thus,  greater  realism  should 
be  introduced  into  the  programs  at  the  cost  of  increased  complexity  and  bookkeeping  problems.  Some 
calculations  utilizing  this  technique  are  compared  with  the  original  Berry  and  Herman  method  in  Figs. 
58  and  59.  We  see  that  for  long  daytime  paths  the  two  progtams  give  similar  predictions,  but  for  pas¬ 
sage  through  the  day-night  terminator  there  are  obvious  differences. 

Because  both  the  multimode  and  multihop  methods  have  complimentary  advantages  when  applied 
to  longwave  propagation  predictions,  both  methods  should  be  applied  to  problems  whenever  either  one 
is  doubtfull 

Future  Needs  for  Coverage  Predictions 

Progress  in  the  accurate  prediction  of  long  wave  radio  fields  is  strongly  tied  to  our  ability  to  moni¬ 
tor  and  model  the  lower  ionospheric  D-region  and  to  accurately  map  ground  conductivity  of  the  earth  as 
a  function  of  frequency.  A  medium-sized  campaign  of  aircraft  flights  and  long  term  ground  based 
receivers  could  reduce  our  uncertainties  in  many  areas. 
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Fig.  58  —  Field  strength  vs  distance  predictions  for  Naval  LF 
Broadcast  Station  Greece  on  59  kHz  for  July,  3.5  kW  radiated 
power,  due  north  at  0000  UT  Solar  zenith  angle  at  midpoints 
(3500  km)  is  62°.  Path  is  quite  homogeneous  daytime,  so  old  and 
new  methods  of  prediction  give  similar  results. 


Fig.  59  —  Field  strength  vs  distance  predictions  for  Naval  LF 
Broadcast  Station  Greece  on  59  kHz  for  July.  3.6  kV*  radiated  power, 
due  north  at  0000  UT.  Solar  zenith  angle  at  midpoint  (3500  km)  is 
87°  or  daytime  condition;  at  2500  km  the  solar  zenith  angle  is  96c  or 
nighttime  condition.  The  old  method  interprets  the  path  as  all 
daylight;  the  new  method  has  reduced  attenuation  for  the  nighttime 
portion  of  the  path  and  appropriate  daytime  attenuation  for  the 
daytime  portion 
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-The  possibility  for  making  a  worldwide  D*region  sensing  and  monitoring  system  is  very  appealing. 
Numerous  systems  can  produce  D-region  ionospheric  information,  and  special  systems  can  be  reudily 
designed  and  implemented  to  serve  this  purpose.  For  example,  the  auroral  photographs  currently  avail¬ 
able  from  the  DMSP  satellite  can  be  correlated  with  actual  propagation  calculations  end  measurements 
to  get  a  model  of  the  aurora-dominated  situation.  Data  from  DMSP  type  sensors  working  in  various 
spectral  bands  could  detect  regions  of  high  electron  precipitation  to  provide  real-time  feedback  to  com¬ 
munications  traffic  managers  to  permit  efficient  mode,  transmitter  location  and  frequency  selection. 
Numerous  inexpensive  ground  systems  can  be  implemented  to  provide  auxiliary  data.  A  prime  exam¬ 
ple  of  a  ground  based  D-region  ionosphere  sensor  would  be  the  large  Arecibo  Thomson  Scatter  Radar. 
It  gives  local  electron  density  profile  information  from  60  km  up  to  1000  km.  Performing  a  series  of 
experiments  with  this  facility  tied  in  with  OTH  HF  radar,  partial  reflection  sensors,  riometers,  and 
microwave  and  airglow  monitors  should  greatly  enhance  our  knowledge  of  the  D-region  and  of  its  prop¬ 
erties  relevant  to  long  wavelength  propagation  prediction.  A  second  series  of  experiments  are  the 
auroral  zone  using  the  University  of  Tromso’s  Partial  Reflection  Experiments  and  EISCAT  Thompson 
Scatter  Radar  facility  at  Tromso,  Norway,  would  give  additional  validity  to  the  polar  ionospheric  models 
and  help  with  the  validation  of  a  polar  imaging  system.  The  data  from  the  HILAT  satellite  discussed  by 
Huffman  et.  al..  [521  are  another  example  of  valuable  satellite  produced  information. 

Dispersion  and  Multipath 

Because  wide  bandwidth  or  short  pulse-length  modulation  can  be  used  for  VLF  communication 
systems,  it  is  important  to  consider  how  multimode  propagation  and  dispersion  (end  to  distort  short  sig¬ 
nal  pulses.  References  53  and  54  present  various  methods  for  treating  transient  signals  in  a  dispersive 
channel 

Reference  55  considers  VLF  pulse  stretching  from  the  point  of  view  of  the  waveguide- mode 
theory  of  VLF  propagation  to  show  the  predictions  that  this  theory  makes  about  short-pulse  transmis¬ 
sion. 


When  short  pulses  are  propagated  in  a  waveguide,  the  received  pulse  shape  is  altered  by  both 
multimode  (or  equivalent  multipath)  propagation  and  bv  dispersion  of  the  pulses  propagating  by  indi¬ 
vidual  modes.  Multimode  distortion  arises  when  a  pulse  leaves  a  transmitting  antenna  and  travels  down 
the  waveguide  in  several  modes.  The  different  modes  have  different  phase  and  group  velocities,  so  the 
pulses  arrive  at  the  receiver  displaced  from  one  another  in  time.  Depending  on  how  much  displace¬ 
ment  has  occurred,  the  pulse  can  look  like  one  long  pulse  or  like  a  series  of  two  or  more  pulses.  Such 
multimode  distortion  is  identical,  in  effect,  to  the  classical  multipath  distortion  that  always  arises  when 
a  signal  can  take  more  than  one  path  to  a  receiver. 

Dispersive  distortion  is  independent  of  the  multimode  problem.  A  short  pulse  is  composed  of  a 
spread  of  frequencies  around  a  carrier  frequency  and  thus  occupies  a  certain  bandwidth  A /  that  is 
roughly  given  by  the  reciprocal  of  the  length  of  the  pulse  t0;  i.c.,  A /  —  1/t0.  Since  the  group  velocities 
of  the  components  on  different  frequencies  are  not  the  same,  the  components  arrive  at  the  receiver  dis¬ 
placed  in  time.  This  leads  to  pulse  smearing  due  to  propagation  even  on  a  single  mode. 

It  is  important  to  remember  that  multimode  distortion  and  dispersive  distortion  are  distinct  and 
independent  phenomena  Thus,  propagation  channels  that  have  a  great  deal  of  multimode  distortion 
but  little  dispersion  and  vice  versa,  could  exist. 

The  equations  that  describe  pulse  stretching  due  to  multimode  propagation  in  the  waveguide  are 
straightforward.  The  index  of  refraction  for  the  nth  mode  nn  is  given  by  the  ratio  of  the  speed  of  light 
to  the  phase  velocity  (we  use  the  letter  n  both  for  the  index  of  refraction  and  as  a  mode  index). 

n„  =  clvph  „  (33) 
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The  group  velocity  v,  „  of  a  signal  on  the  nth  mode  is  obtained  from  the  index  of  refraction  and 
its  derivative  as  follows: 

v*,n/c  =  ( n„  +/  dnjdf)~x .  (34) 


The  time  delay  kTm„  between  pulses  arriving  on  different  modes  (the  mth  and  nth  modes)  is  simply 
given  by: 


A  7_ 


a© 


(35) 


Here  a©  ( =d )  is  the  distance  between  the  transmitter  and  receiver.  The  equations  that  describe  the 
dispersion  spread  are  also  obtained  very  straightforwardly  (56).  The  difference  in  arrival  times  Ard,sp 
between  the  components  on  the  different  frequencies  on  a  single  mode  are  calculated  for  a  frequency 
difference  A/: 


Ard„p  =  fl0 
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v,</> 
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v#  (/  +  A/) 


(36) 


Collecting  terms  and  introducing  the  initial  pulse  length  t0  =  1/A/,  we  obtain 

a©  dv. 


At  =  */ 

a  'diip  2 

v/ 


Vf  (/+A/)-v,(/) 


A  / 


Tov; 


df 


The  fractional  dispersive  spreading  is  given  by 

a©  d 

Tu  rj  V*J  df 


(37) 


(38) 


There  is  a  further  type  of  dispersive  distortion  that  cannot  be  treated  by  these  simple  concepts. 
This  type  of  spreading  occurs  because  of  the  frequency  variation  of  the  excitation  facmrs  and  of  the 
attenuation  coefficients  on  the  various  modes.  The  amplitudes  of  different  frequency  components  of  a 
pulse  are  excited  differently  and  attenuated  differently.  A  pulse  that  was  initially  composed  of  frequen¬ 
cies  over  an  initial  effective  bandwidth  will  appear  to  be  shifted  in  frequency  or  to  occupy  a  smaller 
bandwidth  because  of  the  greater  absorption  of  the  frequency  components  at  one  side  or  the  other  of 
the  occupied  band.  This  type  of  distortion  could  be  very  significant,  but  it  appears  that  it  can  only  be 
treated  by  considering  various  pulse  shapes  and  performing  Fourier  transforms.  The  procedure  for 
numerically  studying  this  type  of  dispersion  is  well  known.  A  pulse  shape  v,  (r)  at  the  transmitting 
antenna  is  given  by  some  modulation  function  A < i)  times  a  basic  sine  wave  exp  /u>0'  on  the  frequency 

<a0: 

v, (/)  =  A  (;)  exp  /o>„r .  (39) 


This  pulse  can  be  Fourier  analyzed: 

1  /'+'" 

v, (w)  = - -  A  (/)  exp  /oinf  rexp  -iwi  <it 

( 2  7t  )  ' 2  J 


(40) 


The  received  pulse  v,  „(oj)  on  the  nth  mode  will  be  given  by: 

e~ikaH 

v,  „(<o)  =  v,(cu)  - 

a  (©sin©) 


i  ?  n 


(41) 
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where  Vn  is  the  nth  term  in  the  summation  of  (2);  i.e., 

y„  -  sin3/20„  A„  exp  {iw/4+ik0  d  sin  0„) .  (42) 


The  received  pulse  must  then  be  transformed  back  to  the  time  domain  to  obtain  the  pulse  shape  as  a 
function  of  time: 


(r)  = 


1 

(2tt)i/2 


(43) 


If  we  desire  the  total  received  pulse  shape,  we  merely  sum  the  pulses  on  the  various  modes: 

v,(f>-  t  (44) 

n  ■"  1 


We  could  just  as  well  obtain  v,(r)  by  performing  the  summation  over  the  Vn  before  the  integration 
over  a).  We  would  then  have  obtained: 


vr  (  / ) 


— 1 —  c  e«*i 

(2t r),/2 


v.(o.)e'*aH 
a  (0  sin  0)|/2 


doj  . 


(45) 


The  equivalence  of  (44)  and  (45)  indicates  that  it  makes  no  difference  whether  we  consider  the  pulses 
as  carried  by  different  modes  or  whether  we  just  consider  the  existence  of  the  complicated  complex 
transfer  function: 


H(u>)  = 


,-iAaO 


a  (0  sin  0) 


(46) 


It  is  easier  to  think  of  pulses  ttaveling  by  different  modes  or  hops,  so  we  prefer  to  use  this  separate 
mode  approach  rather  than  the  transfer  function  approach. 

Dispersion  Calculation 

The  isotropic  exponential  ionosphere  model  with  a  0=  0.5  and  heights  of  70  and  90  km  may  be 
used  to  represent  daytime  and  nighttime  ionospheric  conditions.  From  data  given  in  Ref.  23,  one  may 
calculate  the  index  of  refraction  (Fig.  60)  assuming  an  infinite  ground  conductivity.  The  derivative  of 
the  index  of  refraction  can  be  determined  numerically  (Fig.  61);  and  thus  the  group  velocity  as  a  func¬ 
tion  of  frequency  is  obtained  (Fig.  62).  Using  these  group  velocities,  the  delays  between  pulses  on 
various  modes  were  calculated  The  multimode  interference  will  normally  have  the  greatest  observable 
effect  when  two  or  more  modes  are  approximately  of  equal  strength.  The  region  of  frequency-range 
space  in  which  two  or  more  modes  are  within  4  dB  of  each  other  was  determined  by  studying  the 
graphs  of  signal  strength  vs  distance  given  by  Ref.  24,  which  are  based  on  the  parameters  given  by  Re! 
23.  The  group  delays,  the  time  of  arrival  between  pulses  on  the  first  and  second  modes,  and  between 
pulses  on  the  second  and  third  modes  were  then  computed.  These  delay  times  are  shown  in  Tables  3 
and  4,  and  two  time  delays  are  given  for  most  frequencies.  The  upper  time  delay  is  between  the  first 
and  second-mode  pulses;  the  lower  time  delay  is  between  the  second-  and  third-mode  pulses.  For  some 
frequencies,  Ref  23  gives  the  phase  velocity  for  only  two  modes;  and  for  these  cases  only  one  delay 
time  is  presented.  Delay  time  values  are  given  only  for  the  region  for  which  two  or  more  modes  ha’  e 
amplitudes  within  4  dB  of  each  other.  Delay  lime  values  at  other  ranges  can  be  easily  calculated  by 
proportionality.  In  the  regions  given,  two  pulses  may  be  separated  by  the  time  delay  given  on  the 
upper  line,  or,  if  the  two  dominant  pulses  are  arranged  so  as  to  cancel  out  to  a  great  extent,  the  pulse 
remnant  of  the  cancellation  of  the  first  two  pulses  would  be  followed  by  a  pulse  traveling  on  another 
mode. 
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Fig  62  —  Group  velocity  as  a  function  of  frequency. 
Solid  line,  nighttime,  dashed  line,  daytime 


Table  3  —  Daytime  multipa’th  time  delays  in  milliseconds 
(/3  —  0  S,h  —  70km) 


Frequency, 

kHz 

Range 

,  Mm 

0.5 

10 

1.5 

2.0 

4.0 

6.0 

16 

0.039 

18 

0.031 

0.061 

20 

0.024 

0  047 

0.071 

0.052 

0.104 

0.156 

22 

0.018 

0  036 

0.054 

0.073 

0.042 

0.084 

0.126 

0.168 

24 

0  014 

© 

© 

0.056 

j 

0.036 

0.107 

0.142 

26 

0.009 

0.038 

0.076 

0.031 

0.123 

0.246 

28 

0.008 

0.066 

0.101 

1 _ 

0.025 

0.202 

0.302 

Sole.  Where  two  values  of  delay  arc  given,  the  upper  value  is  the  delay 
between  pulses  on  the  firsi-  and  second-order  modes,  and  ihc  lower  is  ihe 
delay  belwcen  pulses  on  ihc  second-  and  third-order  modes 
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Table  4  —  Nighttime  multipath  time  delays  in  milliseconds 
(p  -  0.5,  h  -  90  km) 


Note.  Where  two  values  of  delay  arc  given,  ihe  upper  value  is  ihc  delay  between  pulses  on  the  first-  and  second-order 
inodes,  and  the  lower  is  the  dela*  between  pulses  on  the  second-  and  third-order  modes 


The  dispersive  time  delay  fraction  Ar/r0  from  Eq.  (38)  was  also  calculated  for  an  initial  pulse 
length  of  2  ms  received  at  6000  km  from  the  transmitter.  The  results  are  given  in  Figs.  63  and  64  We 
see  that  the  spreading  is  small,  typically  only  a  few  percent. 

Results  of  Dispersion  Calculations 

One  observes  that  the  time  delays  between  the  second-  and  third-order  mode  pulses  are  usually 
greater  than  the  delay  between  the  first-  and  second-order  mode  pulses  During  the  daytime  the 
greatest  delay  between  the  first-  and  second-mode  pulses  in  the  interference  region  is  approximately  0.1 
ms  at  28  kHz  and  6-Mm  range  The  greatest  delay  between  the  second-  and  third-mode  pulse  is 
approximately  0.3  ms  at  the  same  frequency  and  range.  At  nighttime  the  interference  region  is  shifted 
to  lower  frequency.  The  greatest  time  delay  between  pulses  on  the  first  and  second  modes  is  approxi¬ 
mately  0.18  ms  at  a  frequency  of  20  kHz  and  12-Mm  range.  The  time  delay  between  the  second-  and 
third-mode  pulses  is  approximately  0.76  ms  at  the  same  frequency  and  range. 
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Dispersion  Conclusions 

Theoretical  estimates  were  obtained  of  the  severity  of  muitimode  and  dispersive  distortion  in  the 
VLF  propagation  channel.  It  appears  that  multimode  propagation  will  cause  greater  pulse  distortion 
than  dispersion  will.  Even  at  long  distance  at  night,  the  multimode  propagation  introduces  appreciable 
pulse  distortion  on  some  frequencies  for  pulse  lengths  shorter  than  —  1  or  2  ms.  since  the  maximum 
delay  time  is  of  the  order  of  0.7  ms.  During  the  daytime,  however,  the  delays  are  smaller,  only  about 
0.3  ms,  and  the  pulse  stretching  is  also  smaller.  From  the  present  calculation,  it  would  appear  that  a 
communicator  can  choose  an  optimum  frequency  for  transmission  to  avoid  the  regions  of  frequency- 
range  space  given  in  Tables  3  and  4  and  thereby  avoid  both  problems  of  multimode  reception  time 
delays  and  of  modal  interference  nulls. 
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Fig.  63  —  Fractional  lime  delay  dispersion  vs  f'equc'ey  Tor  in.ee  modes 
for  a  distance  D  -  iti  of  6  Mm  for  n’-.htlime  ( i  —  0  f>.  6-90  km) 
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Fig.  64  —  Fractional  time  delay  dispersion  vs  frequency  for  three  modes 
for  a  distance  D  •>  aO  of  6  Mm  for  daytime  (/3  —  0.5,  h  =  70  km) 

MODULATION  TECHNIQUES  FOR  LONGWAVE  SYSTEMS 

Most  of  the  time,  the  longwave  communication  channel  was  dominated  by  on-off  keyed  (OOK) 
systems  transmitting  Morse  code.  Frequency-shift  keyed  (FSK)  systems  were  introduced  to  provide 
data  transfer  at  higher  speed.  The  performance  of  these  systems  is  summarized  in  Ref.  2.  Recently, 
minimum  shift  keying  (MSK)  systems  have  been  introduced  to  provide  still  more  efficient  use  of  the 
limited  frequency  spans  of  each  longwave  communications  channel.  The  following  mathematical 
development  is  intended  to  provide  an  understanding  of  the  newer  MSK  systems. 

Elementary  Signal  Theory 

Figure  65  shows  schematically  the  outlines  of  a  long-wavelength  communication  system.  Figure 
66  shows  a  conceptual  version  of  a  receiver  and  demodulation  portion  of  Fig.  65.  The  operation  of  the 
integrator  portion  of  the  equipment  in  Fig.  66  is  to  sample  the  received  waveforms  h(i)  and  to  multiply 
them  by  the  desired  waveform  and  /;2(f)  and  to  integrate  over  a  time  equal  to  a  bit-keying- 

eltment  period  Tg,  centered  about  the  time  r,  the  midpoint  of  the  /  th  expected  bit  interval.  The 
transmitted  waveform  is  constructed  as  a  linear  superposition  of  h\(t)  and  /; 2 ( z )  according  to: 

hit)  =  X,  h\{t)  +  Y,  h2(l).  (47) 

The  outputs  of  the  integrators  yield  the  values  of  X,  and  T,  that  contain  the  transmitted  information, 

x,  =  (Tb)-'  h\(t)hii)  dt 


(48) 
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NOISE 


VLF  PROPAGATION  CHANNEL 


Fig.  65  —  Model  of  the  communication  system 


Fig.  66  —  Model  of  the  receiver-demodulator 


and 

Y,  -  W  X,-r^2  M'>A('>  dt, 

as  long  as  h\ (t)  and  h2(t)  are  orthonormal  to  each  other: 

( T„ )_1  f'V™  h,(t)  h,U)  dt  -  6,. 

Here  S0  is  the  Kronecker  delta  function,  equal  to  zero  for  /  ^  j  and  equal  to  1  for  /  —  J. 


Correlation  o/ Orthogonal  Waveforms 


(49) 

(50) 


A  set  of  orthogonal  wideband  waveforms  /» , (f )  and  h2(t)  can  be  generated  by  using  a  baseband 
waveform  hH(t)  to  modulate  both  a  cosinusoidal  and  sinusoidal  wave: 

/i|(f)  -  hB(t)  cosa>0/,  (51) 

/i2(/)  -  hB(t)  sin  (D0t,  (52) 

h(t)  —  hB(t)( X/  cos  u/0t  +  Y,  sin  w0t).  (53) 

The  values  of  X ,  and  K,  give  the  strength  of  the  cosinusoidal  and  sinusoidal  components.  A  correlation 
receiver  can  then  be  made  to  extract  the  value  of  X and  Yt\ 

XI j  4"  7 g/  2 

T  hB(t)(Xi  cos  u>0:  +  Y,  sin  u>0t)  hB(t)  cos  (u0t  dt  (54) 

and 

Yi  -  2(Tg)~'  I  hB(l)  (X,  cos  w0t  4-  Y,  sin  u0t)  hB(t)  sin  u>0/  dt.  (55) 

*"/  ‘Bu 

A  propagation  delay  time  of  rp  modifies  the  correlator  outputs  to  give  new  ones  that  we  call 
Xj(Tp)  and  Y/(tp): 

X,(rp)  -  2(TB)~i  J  “  hB(t  +  rp)  [X,  cos  w0(/  +  t p)  +  Y,  sin  w0(t  +  t p))  (56) 

x  hH  (t )  cos  u j(,f  di 
and 

Y,(Tn)  =  2(  Ty)~l  S,’-]  “t  +  TpHlA";  cos  w0(t  +  7 ■  p)  +  Yi  sin  w0(/  +  rp )  ]  (57) 

x  hH{t)  sin  w, ,r  dt. 
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We  find  that  the  vector  R,(rp)  -  (*,(tp),  T,(tp))  is  approximately  the  phase-routed  version  of 
R,  -  ( Xt,Y We  have 

r  ‘ i *  rBn 

X,UP)-  Ti'iX,  cosu.0 rp  +  Y,  sino,0Tp)  J,rV2  M*>  M'  +  dt  (58) 

and 

/■  'i+'a^ 

v2  hB^t)hBU  +  rp)  dt.  (59) 


Thus  we  can  write  this  equation  in  vector  form  if  R/  —  (A/,  Y2)\ 


r'l+TJi 

w  -  '■«-|[XrV, 


hBU)hB(t  +  rp)dn  I 


cos  w0Tp 


sin  tii0Tp 


sin  fa»orp 

COS  «U0 Tp 


R, 


(60) 


The  magnitude  of  the  vector  R/(rp)  “  l A”/ (-rp ) 2  +  Y,(Tp)2Vn  is  given  by 

|R,(T,)|  -  VI*/1  S‘'-TBBn 
The  phase  of  the  vector  R,(rp)  is  given  by 

<t>/ (rp)  -  arg(X,(rp )  4-  iY,(rp))  -  arg(X,  +  iY,)  o>0Tp 


(61) 

(62) 


Thus  the  magnitude  of  the  R,(rp)  vector  varies  slowly  with  rp.  and  the  phase  d>,(rp)  changes 
rapidly. 


Correlation  of  Other  Orthogonal  Waveforms 

A  more  complicated  set  of  orthogonal  waveforms  can  be  generated  by  using  two  orthogonal 
baseband  waveforms  hB](t)  and  hB2(t)  such  that  for  /,  j  -  1,  2 

(rfl)->  hBl(t)hBl(t)  dt  -  8„.  (63) 

4/  ,  gli 

where  f>„  is  the  Kronecker  delta  function  The  transmitted  waveform  6(r)  is  given  by 

h(t)  =  XihB[(t)  cos  w0t  +  y,/»fl2(/)  sin  a,0t.  (64) 

To  extract  the  value  of  X,  and  Y,  from  this  waveform,  we  need  a  correlator  receiver  that  does  the  fol¬ 
lowing  integrations: 

X,  -  2 ( 7*0 > “ 1  U,Wf>cosa,„r+  Y,hB2U)  sin  w0r][/ifl,(/)+/ifl2(/)l  cos  u>0/  dt  (65) 

• /lr'u 11 


and 


Y,  -  2(  Tfl)'1  f  '  T®,2  \X,hBi(t)  cos  wq/  +  Y,hB2(t)  sin  u>0tHh„,lt)+hB2(t)l  sin  unt  dt. 
he  time  of  arrival  is  changed  to  t  =  /  +  rp,  we  have 

jr;<Tp)  -  2 ( )  *  J°  r^2  [X,hBl(t  +  r p)  cos  oiq(i  +  r„)  +  K,/7tf2(r  +  rp)  sin  w0(/  +  t„)I 


(66) 


x  [hB,(t)  +  hB2{t ))  cos  at 0l  dt 


(67) 


and 


Y,( Tp)  -  2(7'fl)_l  X':;7;,  +  Tp )  cos  W0(r  +  tp)  +  Y,hB1it  +  rp)  sin  w0(r  +  tp>) 


x  (6sl(r)  +  sin  tu0/  J/. 


(68) 
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This  gives 

Xi(rp)  -  X,  cos  (Oqt p  Tg  *  f  hB\(t  +  T  p)[hB{(t)  +  hgjit) i  dt 

+  Y sin  <0$t p  J hgj(t  +  Tp)\hg\(t)  +  hg}(t))  dt  (69) 


and 


Y/( T„)  -  -A';  sin  (oqt p  7 i» +  T dt 

+  Y i  cos  «oT/>  Tj  lJt  T*n  haiit  +  t„)  {A$|(r)  +  hBJ(t)  1  dt. 

(70) 

if  the  two  baseband  modulation  functions  are  sufficiently  orthogonal  that 

.  r'i*Ten 

Tb  Si,-  Tf/1  h*,K' 

,)h aj(f)  dt  0 

(71) 

and  sufficiently  similar  that 

v  s;;Z  ^>(/ + vw'>  d>  -  tbi  s, 

/,+  Tpn 

rTg, j  hB2(t  +  tp)  hB j(<)  dt, 

(72) 

we  again  obtain 

/* 

Tb/1 

A 'i(tp)  -  (X,  cos  w0tp  +  Y,  sin  w 0tp)  Tb  J( 

,  -re/2  W'  +  rp)hB |(»)  dt 

(73) 

and 

C  ‘i+  Tgll 

Yi(rp)  -  (-  X i  sin  u0Tp  +  Y,  cos  w0rp)  T£ 1 

Ji,-Ttn  *W/+Tf)A*i(*)  dt- 

(74) 

When  the  more  complicated  correlators  of  Eqs.  (65)  and  (66)  are  used,  and  the  orthogonal  waveforms 
satisfy  the  conditions  of  Eqs.  (71)  and  (72),  the  type  of  waveform  of  Eq.  (64)  gives  a  behavior  similar 
to  that  of  Eq.  (53). 


Modulation  Methods 

The  values  of  R,  can  be  modified  at  the  transmitter  to  initiate  the  communication  of  information. 
In  amplitude  modulation  the  magnitude  of  R/  is  varied,  in  phase-shift  keying  the  phase  angle  <t>i  is 
modulated.  Thus,  in  general,  information  may  be  transmitted  by  changing  the  transmitted  vector  R,  in 
some  prescribed  way.  The  stability  of  the  received  vector  R ,(tp)  in  amplitude  and  phase  is  important 
in  the  performance  of  communication  systems,  in  the  realistic  communication  channel  the  values  of 
R,(r^)  will  also  depend  on  the  noise  present.  In  this  study  we  neglect  noise  and  simply  study  the 
effects  of  propagation  on  R/(t,,). 

MSK  Signals 

In  the  field  of  VLF  communications,  minimum  shift  keying  (MSK)  is  a  commonly  used  type  of 
spread-spectrum  modulation.  MSK  is  mathematically  described  in  terms  of  a  transmitted  waveform 
j(r),  a  baseband  modulation  u(r),  and  a  chip  waveform  p(t),  as  follows: 

s(r)  -  Re  u(r)e'“'°'j,  (75) 

, 

u(t)  -  £  hr,/,  p  (I  -  nil  -  i  Y,  Qn  P  (r-  7~/2  -  nT)  ,  (76) 

’  kil-lfll 

p{l)  “  0,  |  /  |  >  77 2.  '  -  T/2  <  I  ^  7V 2,  (77) 
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This  is  a  specific  version  of  the  signal  described  in  genera!  by  Eq.  (64),  where  the  orthogonal  baseband 
waveforms  hBlit )  and  /ifl2(/)  are  given  by 


hB\(t)  -  y/T 


/!“  0 


/„/>(/  —  nT) 


(78) 


and 


hg  j(r)  =  y/T 


1  QnP^t  ~  T/2  -  nT). 

n— 0 


(79) 


In  Eqs.  (76),  (78),  and  (79)  /„  and  Q„  are  pseudorandomly  chosen  integers  of  value  +1  or  -1. 
It  is  assumed  that  both  the  transmitter  and  receiver  possess  synchronized  pseudorandom-number- 
sequence  generators  that  provide  the  /„  and  Q„  sequences.  The  frequency  «0  is  the  carrier  frequency 
of  the  broadcast.  The  time  interval  T  is  the  duration  of  a  chip.  Since  two  chips  are  transmitted  during 
each  chip  duration,  the  chip  length  rc  is  77 2.  The  MSK  modulation  method  gives  rise  to  a  constant- 
amplitude  waveform  with  continuous  phase  as  long  as  X ,  and  Yt  are  equal  to  ±  1.  The  bandwidth  of 

the  transmitted  signal  is  inversely  related  to  the  chip  period.  The  value  of  the  waveform  si/)  depends 

on  one  value  of  /„  and  one  value  of  Q„  simultaneously  except  at  the  instant  when  a  transition  is  being 
made  from  one  value  of  /„  or  Q„  to  the  next.  At  the  transition  time  the  amplitude  of  the  waveform 
undergoing  transition  is  zero,  because  the  pit)  function  goes  smoothly  to  zero  at  the  beginning  and  end 
of  its  range  At  that  instant  the  value  of  sit)  depends  only  on  the  other  integer  Qh  or  /„.  Considering 
the  autocorrelation  function  Bi t)  for  the  MSK  waveform  sit): 

Bir)  -  lim  fr  sit)sit  +  t)  dt.  (80) 

/'-«  T  Jo 

The  use  of  sit)  from  Eq.  (75)  makes  Bir)  nonzero  as  long  as  M  <  T.  When  Irl  >  T.  all  the 
contributions  to  Bir)  are  from  terms  that  have  the  products  Qn  lm ,  Q„  Q„.  and  /„  /„,  with  n  m. 

These  terms  arc  as  often  4- 1  as  they  are  -1  for  an  orthogonal  or  uncorrelated  set  of  [Qn  1  and  \ln  )  ; 

therefore  a  summation  based  on  these  values  should  remain  small.  A  large  value  for  the  autocorrela¬ 
tion  function  occurs  only  for  It |  <  T,  when  /„J  and  Q „2  terms  contribute.  Since  I,}  =  Q, j  =  I,  these 
terms  are  easily  evaluated.  The  correlation  vector  R/iTn)  is  closely  related  to  the  autocorrelation  func¬ 
tion  Bir).  At  very  low  frequencies  the  electromagnetic  waves  may  be  considered  to  propagate  by  dis¬ 
tinct  waveguide  modes  or  hops.  The  time  delays  between  dominant  modes  and  hops  are  small  — less 
than  1  ms  (55) —producing  the  coherent  overlap  of  one  chip  with  its  time-delayed  image.  At  some 
receiver  locations  the  interference  of  the  two  waves  with  each  other  results  in  "selective"  fading.  Such 
fading  is  called  ‘‘selective”  because  it  occurs  only  for  frequencies  for  which  the  interference  condition 
is  satisfied.  At  other  nearby  receiver  locations  or  frequencies  the  time  delays  will  be  such  as  to  produce 
two  pulses  that  reenforce  each  other.  The  propagation  channel  can  be  simulated  by  a  multipath  model 
in  which  h,it)  is  the  transmitted  wave,  and  the  received  wave  is 

s 

hRit)  =  £  4,  hi  it  —  t,),  (81) 

/-i 

where  r,  is  the  propagation  time  for  the  / th  mode  or  hop  and  A,  is  the  amplitude  for  the  / th  mode  or 
hop.  In  the  case  of  earth-ionosphere  waveguide  propagation,  sometimes  the  waves  have  distinctly 
different  group  and  phase  velocities.  Therefore  we  must  consider  both  the  phase  and  group  delay  times 
t j,1*  and  TjJ1'.  The  phase  delay  of  the  baseband  modulation  is  governed  by  and  the  group  delay  is 
governed  by  t^'V 

The  MSK  Correlation  Receiver 

The  performance  of  a  correlation  receiver  for  MSK  modulation  and  the  coherent  detection  of  sig¬ 
nal  reception  on  a  sinple-path  propagation  channel  are  readily  evaluated.  The  correlation  vector  R;(rp) 
obluined  when  one  uses  MSK  modulation  is  explicitly  calculated  by  substituting  the  expression  for 
hu\it)  from  Eq.  (78)  into  Eqs.  (73)  and  (74).  The  integral  !h( t„),  defined  as 
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Mr,)  -  Ti'  hBi(t  +  T„)hBl(t)  dt 

1 

2  ^  / 1  +  7" o/  2 

-  TV1  £  X-r  /2  ~  +  T/>  “  nT) 

_ n  I  8 


(82) 


jt-0 

/»'-0 


is  readily  evaluated.  The  limits  of  integration  in  this  equation  extend  far  beyond  the  region  of  variable 
/,  for  which  the  integrand  is  nonzero.  We  can  change  the  integration  variable  to 

T  -  t  -  nT.  (83) 

The  limits  of  integration  are  thereby  changed  to  /,  +  TB/2  -  nT  and  r,  -  TB/2  -  nT: 

Mr,)  -  M  I  /„/„■  X-+/C7/  ^  +  TP  -  "T>  dt.  (84) 

n-0  1  8 

/»'«■() 

Since  />(/')  —  0,  for  |/'|  >  TI 2,  we  need  to  integrate  only  between  - 772  <  r  ^  T/2.  We  obtain 

T~> 

,2  r  +  r/2 

Mr„)  -  X  J.r/j  /»(/')p (f ’  +  t,  +  (/I  -  n‘)T)  dt.  (85) 


n  “0 
/?'— 0 


From  this  formulation  we  see  that  the  integral  depends  only  on  the  difference  between  n  and  n' 
and  not  on  the  individual  values  of  n  and  n'  separately.  Let  us  define 

r  772 

Q(r,)  “  /  +  rp  +  kT)  dt. 


Since  k  =  n  -  n',  Eq.  (85)  becomes 


M rp)  =  Ti'  £  lnInCn.„(rp) 


/i— 0 
/i  “0 


(86) 


*  |  .  |  s 
7~  "  +  I"T 


-  T’s 1  I  I  /„/„-*  C*(t,) 


?- 


=  7-*-'  I  I  Inln.kCk(rp)U , 


n-k . 


n- 0 


f- 


where 


£/„-*  -  1,  ifO  <  »  -  *  <  y  -  1, 


,V 


0,  if  /j  --  k  <  0  or  n  —  k  >  —  1. 


We  then  interchange  the  order  of  summation  in  Eq.  (87)  to  get 


f" 


7-' 


MM  =  Tu'  £  C\(rp)  £  IJn-kUn_k 

71—0 


(87) 


(88) 


Because  of  the  orthogonality  of  the  set  {/„)  with  any  sequentially  shifted  copy  of  itself,  we  can 


write 


i-i 
2 


,v 


£  In  ln-k  Uft-k  -  f  5*  .0  • 
n=0  ‘ 


(89) 
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This  assumes  that  the  binary  sequence  {/„}  has  been  generated  by  a  method  that  has  "good”  autocorre¬ 
lation  properties  for  the  pseudo-noise-coding  application.  We  obtain 


MO  - 


nt;' 


C0(tJ. 


(90) 


To  calculate  C o(rp),  we  substitute  into  Eq.  (86) 

Co(t„)  =  J  T/2  p(t')pW  +  rp)  dt.  (91) 

The  value  of  C0(tp)  is  symmetric  with  respect  to  rp  and  can  be  obtained  by  using  elementary  integra¬ 
tion  methods: 


(92) 


Equations  (73)  and  (74)  become 

X,(rp)  -  ( X i  cosw0rp  +  Y,  sin  u>0tp)  NT; 'C0(tp) 


and 

K,(rp)  -  (- X,  sin  <o0rp  +  Yj  cos  w0tp) NT;' C0(tp). 


In  matrix  form  we  have 


R  ,(r,) 


nt;  T1 


COS  (OftTp 

-sin  w0rp 


sin  ui0Tp 
cos  (u0rp 


R,  C0( tp). 


(93) 

(94) 

(95) 


A  receiver  can  be  adaptively  synchronized  to  the  incoming  signal  by  introducing  a  variable  delay 
tr  into  the  stream  of  baseband  keying  elements  to  generate  a  family  of  ideal  baseband  keying  elements 
hB\  and  Mathematically  this  means  that  the  variable  t  +  tr  is  introduced  into  the  expression  for 
hHl  in  Eqs.  (68)  and  (69)  in  place  of  t.  The  selection  of  the  synchronization  time  tr  can  greatly  affect 
the  receiver  performance  when  the  actual  propagation  delay  time  rp  is  unknown  (or  equivalently  the 
clocks  that  generate  the  chip  streams  are  not  in  perfect  initial  synchronism).  Therefore  we  rewrite  Eq. 
(95)  to  exhibit  the  dependence  of  the  system  performance  on  tr  Wc  obtain  the  following  expression 
for  R/  (tp,t,): 


R  i(tp,tr) 


NTr 


cos  UI0TP 
-sin  ojqt p 


sin  wqt„ 
cos  w0rp 


R /  ^o(T/> 


T*>- 


(96) 


The  magnitude  of  the  components  of  R,(tp,tr)  is  proportional  to  the  in-phase  and  quadrature  com¬ 
ponents  of  the  received  radio-frequency  carrier  wave  because  of  the  multipliers  cos  wnTr  and  sin  wnrP 
in  the  matrix.  However,  the  magnitude  of  R,(rp,Tft)  depends  on  the  propagation  delay  time  rp  and  on 
the  receiver  time  tr  by  the  factor  rp  =  tr.  If  the  receiver  is  out  of  synchronization,  the  factor  can  be 
zero  or  can  be  very  small.  If  the  receiver  is  employing  phase-shift  keying,  any  change  in  the  propaga¬ 
tion  time  t p ,  even  after  the  receiver  is  perfectly  synchronized  on  the  chip  (  0  ~  I t(1  -  tr  I  «  T),  will 
induce  a  corresponding  phase  change  in  the  phase  of  the  vector  R,  =  (tp,tw),  that  is.  the  doppler 
phase  variation  of  the  radio-frequency  carrier  induces  an  equal  doppler  shift  on  the  phase  angle  d>H . 
These  types  of  doppler-signal  phase  variations  arc  often  seen  when  a  mobile  transmitter  or  receiver  is 
used.  In  such  a  circumstance  we  cannot  assume  that  rp(r)  will  be  constant  at  each  instant  within  a  bit 
period  to  permit  the  factorization  shown  in  Eqs  (67)  and  (68).  However,  the  other  manipulations  that 
we  performed  on  the  expression  fot  R  are  all  valid  up  to  this  point.  To  obtain  the  resultant  correlation 
vector  R,{tp,th)  if  the  delay  time  is  slowly  changing  within  a  bit  period,  we  must  average  the  instan¬ 
taneous  values  of  R/[tp(/),tpI  over  an  entire  bit  period: 
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R/(r,,rs)-N7’s-1rl  f/_T 


r  'i+  Tan 

cos  w orp(f) 

sin  to0Tp(r) 

^  <rT an 

-Sin  too ( Z ) 

COS  to0T p ( f ) 

R;  C0[tp(j)  -  trJ. 


(97) 


From  this  expression  one  can  approximately  evaluate  the  effect  on  the  R/(rp(r),T^)  vector  for  various 
conditions  of  Tp(t)  and  synchronization  conditions  for  a  single  path  propagation. 

We  have  reviewed  the  elementary  cases  of  a  stationary  and  moving  transmitter  that  has  a  single 
path  of  propagation,  since  this  is  the  simplest  case.  Linear  superposition  holds  for  the  waves  and  corre¬ 
lation  procedures;  thus  a  result  from  a  case  that  has  several  propagation  paths  can  be  obtained  by 
adding  the  terms  from  single-path  cases. 


Two  Propagation  Modes  and  a  Stationary  Transmitter 


When  a  stationary  transmitter  produces  a  signal  at  a  receiver,  the  signal  comes  from  two  propaga¬ 
tion  paths,  there  are  two  propagation  delay  times  rpl  and  tp1  and  two  (real-valued)  amplitudes  of  the 
two  modes  A,  and  A:.  Alternatively,  wc  may  say  that  each  mode  is  described  by  a  complex  amplitude, 
and  the  magnitude  of  each  complex  amplitude  is  the  real  amplitude  and  the  phase  is  related  to  the 
propagation  lime.  The  voltage  impressed  on  the  receiving  antenna  by  the  received  signal  represents  a 
waveform  that  is  the  sum  of  two  time-delayed  replicas  of  the  initial  waveform.  The  subscript  i  denotes 
the  mode  index.  The  result  is  shown  in  the  following  expression  for  the  total  correlation  vector 

RT(r«): 


R/(t«)  =  Ri(tpj,th) 
/"I 


-  (V7-*1.  ■'  i  A, 


COS  to0Tp  , 
-sin  cj0tp  , 


sin  to0 tpJ 

COS  toflT,,  , 


Q(iy,  -  tr). 


(98) 


The  behavior  of  R/(r«)  as  a  function  of  the  receiver  chip  synchronization  time  tr  can  be  quite  dif¬ 
ferent  from  ihe  behavior  when  there  is  only  a  single  path.  A  splitting  of  the  correlation  peak  can  occur, 
and  instead  of  the  correlation  being  best  when  tpJ  -  tr,  it  could  be  bad  at  this  time  Numerical  exam¬ 
ples  given  later  exhibit  this  circumstance. 


LONGWAVE  FIELD-STRENGTH  VARIATIONS  FROM  AN  AIRBORNE  ANTENNA  1571 


Consider  the  fields  generated  by  an  inclined,  elevated  transmitting  antenna  as  the  antenna  moves 
in  a  circular  path  (orbit).  This  antenna  is  an  idealization  of  a  trailing-wire  antenna  carried  by  an  aircraft 
circling  at  a  constant  altitude. 

As  the  antenna  goes  around  a  circular  path,  the  VLF  waveguide  modes  generated  by  the  vertical 
and  horizontal  components  of  this  inclined  antenna  cause  periodic  variations  in  field  strength.  The 
magnitude  of  these  orbit-produced  periodic  variations  is  a  function  of  the  distance  from  tne  transmitter. 
The  sense  of  the  variation  (i.e.,  increasing  or  decreasing  at  a  given  time)  can  be  opposite  at  points  that 
are  relatively  close  together  on  the  same  radial  from  the  transmitter.  Either  one  or  two  maximum 
values  of  received  signal  strength  per  orbit  can  be  observed.  When  the  antenna  undergoes  a  pei iodic 
change  of  vertically  caused  by  a  gradient  in  wind  velocity  at  the  altitudes  through  which  the  antenna 
passes  or  droops,  the  magnitude  of  the  periodic  variation  in  received  field  strength  depends  on  the 
direction  of  the  receiver  from  the  transmitter. 


Theory 

Reference  16  gives  the  equations  with  which  the  electric  and  magnetic  fields  can  be  calculated 
from  a  vertical  and  horizontal  antenna  in  the  waveguide  mode  theory.  First  we  consider  only  expres¬ 
sions  for  the  vertical  electric  field  strength,  that  use  the  earth-ionosphere  waveguide  model  in  which 
there  is  no  conversion  of  horizontal  electric  fields  to  vertical  fields  at  the  ionospheric  level. 
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£7'(r,  0)  -  /  W«o>l/J/  ds Vf  £  5„i  5F„ 

n 

E["(r,  9)  -  /(Mo/ «0)l/2/  rfa*  cos<£  X  (*,)5®-5f, 

fl 

where 

i/? 

G<1(z)GB(21)ABexp(i7r/4  +  i*0£>S„)  (101) 

Here, 

£,*(/•,  0),  £/“(/■,  &)  indicate  the  vertical  electric  field  strength  at  the  point  r,  9  received  from  a  vertical 
or  a  horizontal  electric  dipole  transmitting  antenna,  respectively  (Fig.  67). 

(miAo)  1/2  's  the  impedance  of  free  space 
I  is  the  transmitting  antenna  current 

ds1* ,  ds*  .dicate  the  length  of  the  horizontal  or  vertical  transmitting  antenna 

S„  is  the  sine  of  the  complex  eigenangle  for  the  nth  TM  mode 

i„(:,)  is  the  impedance  of  the  nth  TM  mode  at  the  height  zs  of  the  transmitter 

d>  is  the  angle  between  the  great-circle  path  to  the  receiver  and  the  end-fire  direction 

of  the  horizontal  antenna 

\  is  the  free-space  wavelength  of  the  electromagnetic  wave 

D  is  the  great-circle  distance  between  the  transmitter  and  the  receiver 

a  is  the  radius  of  the  earth 

Gn  (z)  is  the  height-gain  function  of  the  mh  mode 

A„  is  the  excitation  factor  of  the  mh  mode 

k0  =  2n/k 

The  total  field  £,(/-, 9)  from  an  inclined  antenna  is  the  sum  of  the  fields  from  the  horizontal  and  verti¬ 
cal  components  of  the  inclined  antenna  when  the  inclined  antenna’s  length  ds  is  broken  up  into  its  hor¬ 
izontal  and  vertical  components  dshe  and  ds*  \  they  are  inserted  in  Eqs.  (99)  and  (100): 

E,(r,  9)  =  E*(r,  9)  +  E,he(r,  9) .  (102) 

To  simulate  mathematically  the  effect  of  the  motion  of  an  antenna  trailing  behind  a  circling  air¬ 
craft,  we  construct  a  Cartesian  coordinate  system  with  its  origin  on  the  earth's  surface  below  the  center 
of  the  circle.  The  z  axis  is  normal  to  the  earth's  surface  and  the  y  axis  points  in  the  direction  of  the 
receiver.  In  these  (x,  y,  r)  coordinates  we  locate  the  antenna  originally  at  the  point  (/•„  O  r,).  where  rv 
is  the  radius  of  the  orbit  and  z,  is  the  height  of  the  antenna  above  the  g.ound.  The  receiver  is  located 
at  the  point  (O,  d.  Ot,  i.e  ,  on  the  y  axis  on  the  surface  of  the  earth  a  great-circle  distance  d  from  the 
origin.  A  diagram  of  this  arrangement  is  shown  in  Fig  68.  The  antenna  is  a  short  straight-line  element 
centered  at  (r„  0,z,)  and  inclined  at  the  angle  y  from  the  normal.  Therefore,  dshr  “  ds  sin  y  and  ds * 
“  ds  cos  y.  The  horizontal  projection  of  the  antenna  points  along  the  y  axis.  We  define  {  as  the  angle 
between  the  xaxis  and  the  projection  in  the  x-y  plane  of  the  line  from  the  origin  of  the  coordinate  sys¬ 
tem  to  the  center  of  the  antenna  Since  d  >  >  r,  we  use  the  approximation  <t>  ~  which  states  that, 


_ I _ D/a 

(Dkh)Ui  sin  (D/a) 


(99) 

(100) 
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Fig.  68  —  Diagram  defining  various  angles 
Upper  diagram  shows  orbit  from  above. 
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when  the  horizontal  projection  of  the  antenna  is  parallel  to  the  .v  axis,  this  horizontal  projection  points 
directly  at  the  receiver.  There  is  a  small  error  in  this  approximation,  but  in  this  calculation  it  is  at  most 
of  the  order  of  r0!d~ 0.01  radian.  The  great  circle  path  D  from  the  antenna  to  the  receiver  is  calculated 
according  to 

D  =•  (<y2+r02)l/2 ,  (103) 

which  is  exact  for  a  flat  earth  and  is  a  good  approximation  for  the  spherical  earth.  With  these  values  we 
then  calculate  Er(r,G)  at  the  receiver.  Next  we  carry  the  antenna  smoothly  about  a  circle,  therefore  we 
choose  some  incremental  angle  At;  and  move  the  antenna  to  the  new  coordinates  (x,',ys',zs')  from  the 
old  coordinates  (x5 ,ys,z5)  that  were  (r„,0,z ,)  according  to  the  transformation 

x,  -  x,  cos  A  £  -  ys  sin  A  £  (104a) 

yi  ~  xs  sin  A  £  +  ys  cos  A  £  (104b) 

z,  =  Zj ,  (104c) 


which  is  a  rotation  of  the  antenna  location  by  an  angle  A£  about  the  z  axis.  The  previous  angles  £  and 
<j>,  were  zero  but  now  are  changed  to  £'  and  <£'  angles  according  to 

<t>'  -  <t>  +  A£  =  £  -  £  +  A£  .  (105) 


A  new  great-circle  distance  £>'  between  the  transmitting  antenna  and  the  receiver  is  calculated  according 


D  =  ((x,)2  +  id  -  >S)J) 1/2 . 


When  we  consider  a  static  inclined  ante.ma  that  is  being  carried  in  a  circle,  the  angle  y  of  the  antenna 
from  the  vertical  is  unchanged.  Because  of  the  changes  in  the  antenna  direction  and  in  the  winds  at 
different  altitudes,  we  must  consider  that  the  antenna  could  undergo  a  periodic  change  in  inclination  as 
it  traverses  a  circle.  We  can  take  this  yo-yo  motion  into  account  by  letting  the  angle  y  oscillate  about 
some  average  angle  y0  us  the  antenna  goes  about  an  orbit  by  setting 

y  =  yo  +  Ay  sin  (£  +  {o)  (107) 


where  £0  is  a  phase  angle  that  relates  the  orbital  variation  of  the  end-fire  angle  <t>  with  the  variation  of  ^ 

the  inclination  angle  y  Static  conditions  exist  when  Ay  “  0.  The  vertical  electric  field  for  the  new  aJ 

antenna  location  is  calculated  according  to  F.qs.  (99)  through  (102).  Then  the  antenna  is  moved 
through  a  further  A£  rotation,  and  the  process  is  repeated.  From  the  results  of  these  calculations,  we  ^ 

determine  the  variations  of  amplitude  and  phase  as  the  antenna  proceeds  about  a  circle. 


Calculations  and  Results 

For  our  sample  calculations,  we  chose  the  mode  parameters  of  Galejs  for  20  kHz  transmissions. 
The  numerical  values  that  we  used  are  given  in  Table  5. 


Table  5  —  Constants  Used  in  the  Calculations 


\  |  =  0.4 

G,(r)  =  1 

Az-  1-2 

G:(z)  =  1 

5;  =  1.002  +  0.000689i 

A.(Zj)  =  — 0.02i 

5  2  *=  0.993  +  0.00 138i 

© 

© 

N 

*1 

z  =  0  km 

r 0=4  !:m  • 

G|  (z, )  =  0.965  +  0.259i 

Zj  »  10  km 

G2(z,)  =  0.95 

Radiated  power  =  1  kw 
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After  each  complete  rotation  we  moved  the  receiver  100  km  farther  from  the  transmitter  over  a  dis¬ 
tance  range  of  500  to  5500  km.  The  topmost  curve  of  Fig.  69  shows  the  vertical  electric  Held  strength 
obtained  by  considering  the  first-  and  second-order  modes  from  a  vertical  antenna  as  a  function  of  dis¬ 
tance.  The  fourth  curve  from  the  top  of  Fig,  69  shows  the  maximum  field  strength  as  a  function  of 
distance  from  an  elevated  horizontal  antenna.  Note  that  the  modal  interference  pattern  for  the  hor¬ 
izontal  antenna  is  quite  different  from  that  of  the  vertical  antenna.  Since  we  consider  an  isotropic 
model  here,  these  results  do  not  include  the  effects  of  polarization  changes  on  reflection  from  the  iono¬ 
sphere. 

The  second  and  third  curves  from  the  top  of  Fig.  69  show  field  strength  as  a  function  of  distance 
for  antennas  inclined  with  y  -  45°  and  y  -  80°,  respectively,  for  Ay  -  0°.  The  error  bars  represent 
the  maximum  and  minimum  values  obtained  as  the  antenna  follows  an  orbit.  The  lower  portion  of  Fig. 
69  shows  the  relative  phase  of  the  fields  from  the  vertical  and  horizontal  antennas. 


Fig.  69  —  Amplitude  and  relative  phase  plots  of 
fields  generated  by  vertical,  inclined,  and 
horizontal  antennas. 

Figure  70  shows  in  more  detail  the  variation  of  amplitude  and  phase  for  an  inclined  antenna  (y  *= 
80°)  at  several  ranges.  Figure  71  shows  the  field  strength  variations  with  distance  lor  the  values  ol 
y0  ™  60°,  Ay  =  11°,  and  £()  =•  0°. 

Discussion  and  Results 

Figure  69  indicates  that  the  amplitude  of  the  orbit-produced,  periodic  variations  of  the  received 
field  strength  is  a  function  of  the  distance  from  the  transmitter.  In  Fig.  70  we  see  a  complete  reversal 
of  the  sense  of  the  oscillation  between  4000  and  4200  km.  At  4100  km  the  field  strength  has  two  rela¬ 
tive  maxima  per  orbit.  This  behavior  seems  very  puzzling,  although  it  has  a  natural  explanation.  Fig¬ 
ure  72  shows  a  complex  phasor  representation  of  the  sum  contained  in  hq  (102).  The  total  vertical 
electric  field  strength  is  the  phasor  sum  of  the  field  strengths  generated  by  the  vertical  and  horizontal 
components  of  the  antenna,  and  £,w  (These  are  the  same  as  £/e  and  respectively.) 
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Fig  70  —  Amplitude  anil  phase 
plots  of  fields  as  a  function  of  the 
lutr.  angle  i,  of  the  antenna 
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Fig  71  —  Amplitude  plot  of  field  generated 
by  an  orbiting,  yo-yoing  antenna 
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Generally,  these  phasors  may  have  any  angle  with  respect  to  one  another;  the  angle  may  vary  with 
distance.  The  E,  phasor  represents  an  omnidirectional  radiation  pattern  that  does  not  change  with  orbit 
(as  long  as  y  is  fixed).  The  £,w  phasor  is  not  omnidirectional,  however,  and  varies  as  the  cos  d>  (see 
expression  100);  therefore,  it  is  modulated  by  a  factor  that  changes  from  +1  to  -1  with  the  orbit  In 
Fig.  72  this  is  represented  by  ±  E,H.  The  resultant  field  strength  is  represented  by  a  phasor,  one  end  of 
which  is  fixed  at  the  origin;  the  other  end  slides  along  the  £,w  phasor.  The  sense  of  the  variation  of  the 
resultant  phasor  is  indicated  by  the  curved  trace  on  the  right-hand  side  of  the  figure.  The  field 
strength,  therefore,  goes  first  through  a  maximum  and  then  to  a  minimum  as  the  antenna  orbits.  The 
middle  diagram  of  Fig.  72  shows  the  ErH  phasor  at  90°  to  the  £/  phasor;  this  relation  causes  two  relative 
maxima  in  field  strength  per  orbit.  When  the  angle  between  phasors  rotates  farther,  as  shown  on  the 
bottom  of  Fig.  72,  the  case  of  a  single  maximum  per  orbit  is  reestablished,  but  now  the  vaiiation  is  in 
the  opposite  sense  to  that  shown  in  the  top  diagram  of  Fig.  72. 


Fig  72  —  The  number  of  maximum  field  strength 
values  per  orbit  may  change  with  the  distance  from 
the  transmitter 


From  Fig.  69  we  see  a  gradual  change  in  the  size  of  the  periodic  variations  with  distance.  We  also 
see  that  quite  different  modal  interference  patterns  are  set  up  by  horizontal  and  vertical  antennas  (The 
difference  is  caused  by  the  presence  of  the  impedence  ^„(z5)  in  Eq.  (100)  and  by  its  absence  in  Eq. 
(99)).  The  ratio  of  E,"  to  £/  changes  with  distance,  and  the  nulls  of  the  two  interference  patterns 
occur  at  different  ranges.  This  effect  is  diagramed  in  Fig.  73.  The  lower  poriion  shows  two  horizon¬ 
tally  generated  phasors  from  two  modes  that  combine  to  give  a  resultant  phasor.  Two  different  E" 
phasors  are  shown  to  illustrate  the  relative  rotation  of  the  horizontally  generated  phasors  of  the  first 
and  second  modes  as  the  distance  from  the  transmitter  changes.  In  the  upper  poriion,  two  vertically 
generated  phasors  also  combine  to  give  a  resultant.  Here  two  different  E\  phasors  are  shown  to  illus¬ 
trate  the  relative  rotation  of  the  vertically  generated  phasors  of  the  first  and  second  modes  as  the  dis¬ 
tance  from  the  transmitter  changes.  The  vertically  generated  resultant  is  then  added  to  the  horizontally 
generated  resultant.  The  horizontally  generated  resultant  is  modulated  by  a  cos  <A  term;  therefore  the 
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Fig.  73  —  Changes  in  modal  interference  with  range 
can  change  the  received  amplitude  pattern 


final  resultant  phasor  may  have  one  or  two  maxima  per  orbit,  as  we  mentioned  previously.  As  illus- 
tiated  in  the  upper  diagram  of  Fig  73,  the  relative  rotation  of  the  mode  phasors  can  cause  the  condition 
to  change  from  a  single  maximum  per  orbit  at  one  distance  to  two  relative  maxima  per  orbit  at  another 
distance.  With  the  help  of  diagrams  shown  in  Figs.  72  and  73,  the  behavior  shown  is  Figs.  69  and  70 
can  be  understood. 

When  the  angle  y  changes  with  the  orbit  angle,  as  in  Eq.  (107),  we  find  an  additional  periodic  variation 
that  enhances  or  reduces  the  received  field  strength  variation  caused  by  purely  orbital  motion.  This  is 
illustrated  in  Fig.  71,  where  a  large,  abrupt  change  in  the  field  strength  variation  is  predicted  between 
4000  and  4200  km  and  between  2400  and  2500  km.  For  these  cases,  the  purely  orbital  variation 
changes  its  sense  of  direction,  whereas  the  yo-yo  variation  does  not  change  appreciable  in  size  or  sense. 
On  the  near  side  of  the  null  at  2400  km,  the  two  fluctuations  are  out  of  phase  and  tend  to  cancel  each 
other  (1.3  dB  oscillation).  On  the  far  side  of  the  null  at  2500  km,  however,  the  fluctuations  are  in 
phase  and  reinforce  to  give  a  12  dB  oscillation.  If  we  had  chosen  4o  10  be  180°  instead  of  0°.  the  situa¬ 
tion  would  be  reversed.  Thus  the  variations  are  dependent  on  bearing  as  well  as  on  range.  The  refer¬ 
ence  direction  for  determining  the  bearing  dependence  is  provided  by  the  wind-shear  direction 

For  the  cases  considered,  we  have  assumed  an  homogeneous  isotropic  ionosphere,  i.e.,  no  earth's 
magnetic-field  effects.  For  this  isotropic  condition,  the  vertical  electric  field  E,  comes  from  the 
transverse  magnetic  cTM)  modes.  The  transverse  electric  (TE)  modes  from  the  horizontal  component 
are  generated  with  a  sin  <t>  variation,  but  they  do  not  carry  any  E,  field.  When  the  earth's  magnetic 
field  is  taken  into  account,  however,  we  find  that  the  TE  modes  become  quasi-TE  and  have  an  appreci¬ 
able  E,  field  that  is  due  to  polarization  conversion  at  the  ionosphere.  This  will  create  an  additional 
complication,  that  is  illustrated  in  Fig.  74.  Here,  the  E'e  phasor  is  fixed  in  magnitude,  but  the  £/"’ 
(TM)  and  E?c  (TE)  phasors  are  modulated  by  cos  and  sin  d>,  respectively,  as  the  antenna  goes  around 
an  orbit.  These  three  phasors  must  be  added  to  give  the  resultant  observed  field: 

Er  =  £.  ^(TM)  +  £,'"’(TE)  .  (108) 
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Fig.  74  —  The  additional  complication  that  the  quasi-TE  modes 
may  make  in  the  received  field  variation 

Actually,  there  will  be  E?  (TE)  and  £,,'f  (TM)  Helds,  but  they  are  combined  to  form  a  constant 
E,"  Held.  The  £/*  (TE)  phasor  has  a  sinusoidal  variation.  The  variation  is  shifted  by  just  90°  in  <t> 
angle  from  the  £,v  (TM)  phasor.  Instead  of  the  end  of  the  resultant  phasor  £,  following  the  line  traced 
out  by  the  £/"  (TM)  phasor,  it  will  now  follow  around  an  ellipse  as  4>  goes  from  0°  to  360°,  as  shown  in 
Fig.  74.  As  £*r  (TM)  goes  through  zero,  £/*  (TE)  reaches  its  maximum  and  vice  versa.  It  is  possible, 
of  course,  for  the  ellipse  to  become  either  a  straight  line  or  a  circle,  depending  on  the  relative  magni¬ 
tudes  and  directions  of  £/*  (TM)  and  £,v  (TE)  on  the  phasor  plane. 

For  an  orbiting  horizontal  antenna,  the  homogeneous,  isotropic  ionosphere  model  predicts  that 
the  Held  strength  will  be  zero  for  <t>  -  90°.  However,  the  introduction  of  the  earth’s  magnetic  Held 
changes  this  prediction;  the  contribution  from  the  quasi-TE  modes  has  a  maximum  lit  <£  -  90°.  Since 
the  quasi-TE  and  quasi-TM  modes  generate  different  modal  interference  patterns,  there  should  exist 
some  combination  of  range,  frequency,  and  bearing  at  which  the  ellipse  of  Fig.  74  is  a  circle.  In  this 
case,  an  horizontal  orbiting  antenna  would  cause  no  amplitude  variation  at  the  receiver.  This  special 
case  is  not  likely  to  be  experienced  often,  but  this  possibility  shows  the  diversity  of  the  results  that  may 
be  obtained  from  an  orbiting  antenna. 

We  see  that  the  vacation  of  the  signal  received  from  an  orbiting  antenna  is  quite  complicated, 
even  though  the  antenna  itself  is  a  simple  short  dipole.  The  magnitude  of  the  variation  is  a  function  of 
the  range  and  bearing  of  the  receiver  Large  changes  in  the  size  of  the  variations  can  occur  with  small 
changes  in  range. 

ANTENNA  MOTION  EFFECTS  ON  MODULATION  SYSTEMS  |58| 

With  the  previous  development  of  the  mathematics  it  is  easy  to  show  the  effects  of  rotation  by  an 
inclined  elevated  antenna.  The  horizontal  component  of  the  inclined  antenna  generates  TE  and  TM 
modes,  the  horizontally  generated  fields  have  phases  and  amplitudes  different  from  those  generated  by 
the  vertical  antenna  components.  The  waveguide  modal  interference  pattern  for  an  elevated  horizontal 
antenna  differs  from  that  of  a  vertical  antenna.  This  effect  is  illustrated  in  Fig  69  It  shows  the  fields 
generated  by  various  inclined  antennas  at  20  kHz,  The  amplitudes  of  the  horizontally  generated  fields 
vary  as  the  cosine  of  the  angle  between  the  horizontal  antenna  component  and  the  direction  of  propaga¬ 
tion 


In  simulating  the  effect  of  this  variation  on  the  receiver  output  vector,  we  use  Eq.  (98)  for  the 
received  correlation  vectors,  we  also  use  four  amplitudes  1-4,!,  four  phase  time  delays  it /I,  and  four 
group  time  delays  |rf).  Two  of  the  amplitudes  represent  two  vertically  generated  mode  fields  and  vary 
as  sin  y.  The  other  two  represent  horizontally  generated  TM  mode  fields  and  vary  as  the  product  cos 
£  cosy.  The  resultant  expression  for  '.he  correlator  output  is  the  same  as  Eq  (98)  except  that  in  this 
case  the  correlation  output  is  the  summation  of  the  outputs  for  the  four  individual  modes. 
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Denoting  the  X^Htr)  from  the  vertically  electrically  generated  modes  as  A'/'lvc(rft)  and  the 
A'/,)  (r/j)  from  the  horizontally  electrically  generated  modes  as  .Y/'’ v(rR ),  we  can  write 

-  cos  y  £  ^(')v'  (t*)  +  cos  £  sin  y  £  A'/'1*'  (t„).  (109a) 

i-i  i-i 


Parallel  definitions  can  be  made  for  Y,t(tr ),  Y,{i)',e(rR)  and  >',(',*<’(tr),  and  we  obtain 

Y,t(.tr)  —  c os  y  £  K//)v,!(r«)  +  cos  {  sin  y  X  (r^ ).  (109b) 

/-i  /-i 

The  vector  R7  (rR)  is  a  sum  of  two  vectors  R /"  (rs)  and  R"  (rfi): 

R,W  -  cos  y  R,v'  ( T R )  +  cos  £  sin  y  R l"(rR).  (110) 

The  vectors  R,vf  (rR)  and  R[*  (tr)  can  be  parallel,  perpendicular,  or  opposed  to  each  other,  depending 
on  circumstances  and  the  value  of  tr.  The  resultant  R/  (tr)  vector  can  change  quite  dramatically  dur¬ 
ing  a  complete  orbit.  For  example,  if  R;'  (r R)  and  R/*  (tr)  are  of  the  same  magnitude  and  phase, 
there  would  be  one  large  oscillation  of  |R7(tr)|  per  orb>t  (provided  y  ”  45T  If  R,’"’(rw)  and  R/V(t r) 
are  perpendicular  to  each  other,  there  would  be  two  oscillations  of  |R7(tr)I  with  one  cycle  £.  There  is 
a  mathematical  isomorphism  between  the  phasor  Er  (D)  and  the  vector  R ,r(rR),  since  the  composition 
of  various  phasors  that  form  the  complex  phasor  sum  E,(D)  is  isomorphic  to  the  composition  of  vari¬ 
ous  correlation  vectors  that  form  the  two-dimensional  correlation  vector  sum  R/(tr).  When  both  TE 
and  TM  waveguide  mode  propagation  from  an  airborne  transmitting  facility  are  taken  into  account,  we 
obtain  the  following  relation. 

R/(t  r)  -  R/'Tt*)  cos  y  +  sinylcos  £R,*'(  TM,tr  )  +  sin  £R  I"(TE,tr)\.  (Ill) 

The  total  vector  R ,r(rR)  is  the  sum  of  independently  propagated  vertical  endfire-  and  broadside¬ 
generated  components.  These  components  are  sums  of  components  that  are  propagated  by  different 
waveguide  modes. 

Diagrammatic  Representations 

It  is  instructive  to  view  diagrams  that  illustrate  the  behavior  of  the  various  R;  vectors  as  r R  and  £ 
vary.  The  following  example  illustrates  what  can  happen  in  a  particularly  severe  case  of  modal  interfer¬ 
ence.  Figure  75  depicts  the  variation  of  a  single-mode  correlation  output  R/i,v<'(tr)  from  a  stationary 
vertical  electric  transmitter.  Figure  76  shows  an  output  function  R/7Iv<,(tr)  for  a  signal  coming  slightly 
later  on  a  second-order  mode  and  180°  out  of  phase  with  the  first.  Figure  77  gives  the  resultant 
r;p(Tr),  that  is  the  sum  of  the  correlation  outputs  of  the  first  and  second  modes.  Figure  78  shows  the 
phase  of  the  vector  R'c(tr )  as  a  function  of  tr. 

Figure  79  shows  the  results  R/(rR)  —  cosy  R (vr(rR)  +  sin  y  cos  £R,',',(tr ) ,  where  £  is^ the  angle 
between  the  direction  of  propagation  and  the  endfire  direction  of  the  antenna.  The  angle  y  is  the  incli¬ 
nation  angle  of  the  antenna;  we  take  y  -  45°  and  £  =  0  for  this  example.  Figure  80  shows  the  same 
results  except  for  the  condition  £  *=  180°.  The  conditions  for  £  “  90°  and  £  =  270  are  both  shown  by 
the  pattern  of  Fig.  77,  since  R ,r(rR)  is  composed  solely  of  the  vertically  generated  modes  in  these 
cases.  These  figures  give  a  picture  of  R,7(tr)  at  the  four  cardinal  points  of  an  orbit. 

In  Fig.  80  four  sample  synchronization  times,  designated  rR),  rR2 •  T«j-  and  rR4.  represent  possi¬ 
ble  correlation  limes  that  might  occur  within  a  receiver.  We  look  at  the  variation  of  the  correlation 
vecior  magnitude  |R,7 <-r )  I  and  phase  6/(tRi)  during  a  complete  orbit  for  these  times  Figures  81  and 
82  show  the  phase  and  magnitude  of  R/(tr,).  The  phase  is  constant  in  this  case,  because  R / r ( t w , )  is 
greater  than  and  parallel  to  R /”( r R | ) .  The  magnitude  varies  in  a  mild  way.  In  Figs  83  and  84  showing 
the  phase  and  amplitude  for  the  time  rR2,  the  phase  makes  two  quick  successive  180°  changes  in  the 
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75  —  R/"'f(rs)  vs  tk  for  a  stationary 
vertical  electric-dipole  transmitter 


Fig  76  —  R,(II,,(r)  vs  t,  for  a  stationary  vertical  electric 
dipole  transmitter.  The  phase  of  this  wave  is  180°  opposite 
to  that  illustrated  in  Fig.  75.  The  maximum  value  of 
I  R/j,-(t«)|  occurs  at  slightly  greater  t*  values  than  does 
the  maximum  value  of  I  RI"”(.tr)\  shown  in  Fig.  75. 


Fig.  77  -  R/'(r„>  -  RrMr*)  +  R/'^Mt*) 
shows  an  abrupt  change  of  phase  at  t R  —  0 
and  a  large  shift  in  the  value  of  r*  for  which 
!  R,'(r,  )  I  is  greatest 
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Fig.  81  —  ( 7 k\)  as  the  angle  of  turn  £  of  an  inclined 
antenna  changes  from  -180°  to  +180°  for  the  synchronization 
time  rx,.  This  phase  is  constant  because  the  contribution  of 
the  vertical  antenna  component  always  exceeds  that  of  the 
horizontal  component  for  this  synchronization  lime. 


Fig.  82  —  I  R,r(T*|)  I  as  £ 

varies  from  -180°  to  +180°. 
The  relatively  gentle  variation 
of  I  R/(r*|)  I  occurs  because 
I  Rr<T*|)  |  >  !  R/v (r * ,)  I. 


Fig  83  —  <£,'(t hi)  vs  orbit  angle  £.  For  this 
lock-on  time  there  are  two  rapid  and  large  phase 
changes  near  £  -  180° 
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Fig.  84  -  I  R,r(T#i)  I  vs  orbit 
angle  £,  showing  how  this  function 
goes  to  zero  when  <*>/( tk:)  is 
making  the  abrupt  phase  transition 


neighborhood  of  £  “  1 80°,  and  the  amplitude  goes  through  two  successive  nulls.  Figures  85  and  86 
give  R,r(0).  This  vector  varies  cosinusoidally,  because  there  is  no  contribution  from  the  vertical  elec¬ 
tric  antenna  component.  Figures  87  and  88  show  the  behavior  of  R/r(r/uK  which  is  similar  to  the 
behavior  of  R/^t^j)  except  that  the  nulls  occur  near  £  -='0°  rather  than  £  —  180°.  Figures  89  and  90 
show  the  behavior  of  R/(t«4),  which  is  similar  to  R /r( r ^ ( )  except  for  a  shift  in  orbit  angle  £  of  180° 


Fig.  85  —  4>it( th  =  0)  vs  angle  £.  There  is 
essentially  no  contribution  from  the  vertical 
component  for  this  value  of  tr,  so  the 
phase  is  controlled  by  the  horizontal 
component  of  the  transmitting  antenna 


Fig  86  —  I  K/(r#  =  0 ; j  vs  orbit  angle  £  The 
cosinusoidal  variation  is  caused  by  the  cos  { 
modulation  oT  i he  Horizontal  component  The 
vertical  component  contributes  essentially  zero 
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Fig.  90  —  I  R/(t/*4)  I  vs  ''•■bit  angle  {.  For  this 
receiver  time  the  maximum  correlation  vector 
occurs  for  5  “  180°,  whereas  for  the 

maximum  correlation-vector  magnitude 
occurred  at  r  —  0°. 

Finally,  Fig.  91  shows  an  estimate  of  t Ropl ,  the  optimum  value  of  receiver  correlation  time  for  dif¬ 
ferent  values  of  orbit  angle  £.  This  graph  shows  that  r/n  and  are  preferable  when  !R,(t*|)|  and 
|R,(t*4)|  are  largest.  The  phase  is  constant,  and  the  amplitude  varies  gently.  Relatively  mild 
amplitude  behavior  occurs  for  this  correlation  time  because  I R /v*’(r « , ) I  >  |R/A,(rR|)|.  The  phase  is 
constant  because  R /v<’(r , )  is  parallel  to  R Although  Rthr(rR])  can  change  the  magnitude  of 
R/fj-fl]),  it  does  not  change  the  phase. 


Fig.  91  -  Values  of  that  should 
be  chosen  for  various  orbit  angles  £ 
for  the  case  illustrated  in  Figs.  75 
through  90 

In  the  pieceding  example  there  was  minimum  correlation  near  the  middle  of  the  correlation 
envelope,  and  the  correlation  maximum  was  greatly  displaced  from  the  maximum  point  of  the  correla- 
lion  of  either  individual  multipath  signal.  Thus  the  rather  small  difference  in  arrival  times  r,,  -  t „2 
changes  the  optimum  chip  synchronization  time  r H  by  a  major  amount  This  effect  is  fairly  easy  to 
explain.  The  condition  described  by  our  equations  is  a  severe  case  of  multipath  interference.  Two  sig¬ 
nals  come  to  the  receiver  antenna  offset  in  time  by  rpl  -  -p}.  The  first  and  second  signals  travel 
independently  of  each  other,  so  when  the  first  s>gna!  pulse  arrives,  the  signal  begins  to  build  up  Then 
the  second  signal  arrives  out  of  phase  and  cancels  ou;  the  first  signal  for  a  long  interval  of  time  while 
both  are  strong  Then  the  first  signal  goes  to  zero  because  its  waveform  has  passed  by,  while  the 
second  signal  remains  to  give  an  excitation  of  the  antenna  again.  Thus  the  strongest  portions  of  each 
signaling  element  in  this  case  of  interference  occur  at  the  beginning  and  the  end  of  each  signaling  ele¬ 
ment.  In  the  middie  of  the  element,  the  two  signals  arriving  from  different  paths  interfere,  and  a 
icduction  or  the  field  present  appears. 
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When  a  receiver  tries  to  synchronize  with  the  signal  at  its  strongest  point,  it  must  search  for  syn¬ 
chronization  either  at  the  beginning  of  the  chip  or  at  the  end  of  it  but  not  in  the  middle.  Thus,  if  a  sys¬ 
tem  were  presynchronized  with  an  optimum  r R  for  either  of  the  individual  mode  delays  or  for  an  aver¬ 
age  of  the  two  mode  delays,  it  would  actually  be  poorly  synchronized  for  optimum  operation!  For 
better  operation  three  correlators  might  be  used— one  adjusted  at  the  central  expected  time  of  arrival  of 
the  wave,  and  one  at  each  end  of  the  chip.  If  the  two  waveforms  reinforce  each  other  to  give  the  max¬ 
imum  in  the  middle,  the  central  correlation  peak  will  be  large  and  the  two  side  peaks  will  be  small.  If, 
however,  an  interference  condition  is  present,  the  middle  correlator  will  be  weak  and  the  two  side 
correlators  will  be  strong. 

Numerical  Calculations 

Procedures 

Calculations  will  follow  of  some  values  of  R,r  for  cases  that  are  inspired  by  realistic  mode-theory 
calculations  of  the  phase  and  amplitude  of  VLF  waves.  For  the  calculated  fields  shown  in  Fig.  6,  the 
waveguide  mode  constants  are  given  in  Table  5. 

These  constants  were  used  in  the  following  equations  to  generate  the  field-strength  values  shown 
in  Fig.  69: 

Er(r.e)  =  /W«o)'/2/  ds”  £  SyFn  (112) 

n 

and 

£,v(r,0)  =  /Wt),/J/  dshf  cos  <t>  n(zs)Sn0SFn,  (113) 


where 


(Da/i)1'2 


D/a 

sin  (D/a) 


G„(z)Gn(z,)A„  exp  (in/ 4  +  ik0  DS„). 


These  equations  are  simplified  versions,  since  they  contain  only  TM-mode  constants  of  an  isotropic' 
ionosphere.  We  calculate  R/',V1'  for  each  mode  according  to 

.V/,Kf  (t*.{)  -  jV|/4,v<’|  cos[w0t /’■'"(£)]  C  (r  *•«(*)  -  r«)  (115) 

and 

T/"'f  (t*,£)  =  jV|/<,v<,|  sinlw0T /’■*'(()]  C  (r  *•'"(£)  -  t«),  (116) 

where 

/V  is  a  normalizing  constant, 

A,vr  is  the  amplitude  of  the  ah  TM  mode  generated  by  the  vertical  component 

of  the  transmitting  antenna, 

o»o  is  the  angular  frequency  of  the  wave, 

t /*■“•({)  is  the  delay  time  of  the  phase  of  the  wave  when  the  aircraft  is  at  angle  £  of 
its  orbit  foi  the  /th  waveguide  mode,  and 

Tf'r(c)  is  the  group  delay  time  of  the  wave  when  the  aircraft  is  at  angle  c  of  its 
orbit  for  the  /th  waveguide  mode. 
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The  amplitude  is  given  by 

A,vr  —  5,'  5G,(r)  G,(ztU,  exp (in/r  +  i  k0  DS,).  (117) 

The  phase  delay  time  r /’■'"(£)  is  obtained  by  evaluating  the  phase  of  A,: 

w0  -  arg(/j,Vf).  (118) 

The  group  delay  time  r*Vf(£)  is  obtained  by  evaluating  the  derivative  of  the  phase  of  A  \ 

T*  V,({)  “  ^[arg(/<'V,)l-  (ll9) 

The  path  length  D  from  the  transmitting  antenna  to  the  receiver  location  is  a  function  of  the  orbit 
angle  £  and  the  orbit  radius  r  according  to  the  approximation 

D  =  D0  +  r  sin  £,  (120) 

where  D0  is  the  distance  from  the  receiver  to  the  center  of  the  orbit.  To  calculate  R/v'(4,t#).  we  per¬ 
formed  summations  as  indicated  in  Eqs  (108)  through  (110). 

We  calculate  the  horizontally  electrically  generated  fields  in  an  analogous  manner: 

X,U)he  (t„,£)  =  N\A*\  cos  [cu0t,  *(£)]  C0  (r «•*({)  -  rR)  (121) 

and 

Y^he  (r r  ,£)  =  N\A* I  sin  U0rf  '*f(£)l  C0  (,■*•*({)  -  t„).  (122) 

The  amplitude  of  the  ,;ih  mode  is  given  as  follows: 

A*  -  Gi(z)G,{zt) A,  exp  (in/ 4  +  ik0  D  S,)  S, 5  A,(zs)  -  A,'e  S,  A,(:u).  (123) 

The  phase  and  group  delay  times  are  given  in  the  same  manner  at  for  the  vertically  generated  modes: 

w0  «  arg(^^)  (124) 

and 

t,4'/"'(£)  =  -^-largM'"')]  U25) 

O(o 

Tables  6  and  7  give  the  numerical  values  used  in  making  me  correlation  calculations. 

Results 

Figures  92  through  97  show  Rf(rR)  and  (1>'(tr)  contours  for  y  values  of  90°.  3  V?.  80',  60".  45°, 
and  0°  with  the  propagation  parameters  shown  to  Table  6.  Figure  98  shows  R"”"Tr/{  t.  the  correlation 
vector  of  the  first-order  mode.  Figures  99  through  104  show  the  correlation  vector  lor  propagation 
parameters  of  Table  7  for  y  value;  of  9u\  85c.  80’,  60’,  48'.  and  0C  The  ''average"  chip  duration  lor 
each  of  these  calculations  was  r,  =  1  ms  The  total  length  ut  each  chip  is  T  =  2r,  =2  ins  The  radius 
of  the  orbit  is  3  km;  the  distance  from  the  center  of  the  orbit  to  the  receiver  is  4000  km. 

Discussion 

Figure  92  shows  the  correlator  performance  for  the  wase  of  a  completely  horizontal  antenna.  For 
this  case  !R,(rK)l  goes  to  zero  at  £  =  90°  and  2 ,,v'.  and  there  are  sudden  1 80c  phase  flips  at  these 
angles.  The  best  correlation  time  appears  tc  be  at  r  H  --  Id..1)  ms,  which  is  the  middle  of  the  correlation 
peak. 
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Table  6  —  Set  of  Quantities  Used  in  Calculating 
Correlation  Vectors 
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Fig.  93(a)  —  Amplitude  behavior  of  R  1(tr)  for  y  85°, 
calculated  with  the  parameters  of  Table  6 
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Fig,  93(b)  —  Phase  behavior  of  R/  (t«)  for  y  -  85°, 
calculated  with  the  parameters  of  Table  6 
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Fig.  94(a)  —  Amplitude  behavior  of  R  1(tr)  for  y 
calculated  with  the  parameters  of  Table  6 


Fig.  94(b)  —  Phase  behavior  of  R/(t«)  for  y  = 
calculated  with  the  parameters  of  Table  6 
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Fig  94(d)  —  Phase  behavior  corresponding  lo  i hai 
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Fig.  95(a)  -  Amplitude  behavior  of  R,t(t«)  for  y  -  60 
calculated  with  the  parameters  of  Table  6 


Fig.  95(b)  —  Phase  behavior  of  R/  (r«)  for  y  -  60 
calculated  with  the  parameters  of  Table  6 
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Fig  96(d)  —  Phase  behavior  corresponding  to  lhal  in  Fig.  96(h) 
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Fig.  97(a)  —  Amplitude  behavior  of  R/  (t#1  for  >  “  0°, 
calculated  with  the  parameters  of  Tabic.  6 
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f  ig  97(b)  —  Phase  behavior  of  R/(t«)  for  y  —  0". 
calculated  with  '.lie  parameters  of  Table  6 
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Fig  97(d)  —  Phase  behavior  corresponding  lo  lhat  in  Fig.  97(b) 
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Fig.  98(a)  —  Amplitude  behavior  of  R/n  (t«)  for  y  0  , 
calculated  with  the  parameters  of  Table  6 


Fig.  98(b)  —  Phase  behavior  of  R/n  (r«)  (or  y  0  . 
calculated  with  the  parameters  of  Table  6 
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Fig  98(d)  —  Phase  behavior  corresponding  10  that  in  Fig  98(h) 
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Fig  99(b)  —  Phase  behavior  of  R /  ( r  „  I  y  -  90‘ 
calculated  with  the  parameters  of  Table  7 
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Fig  100(a)  -  Amplitude  behavior  of  R/ (r*  )  for  y  =-  85". 
calculated  with  the  parameters  of  Table  7 
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Fig  101(b)  —  Phase  behavior  oF  R/ (r*  )  lory 
caleulaicd  with  ihc  parameters  ol  I  able 
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Fig  104(d)  —  Phase  behavior  corresponding  to  that  in  Fig  104(b) 
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Figure  93  shows  the  correlator  performance  for  the  case  of  an  antenna  with  inclination  y  =  85° 
from  the  vertical.  For  this  case  it  is  difficult  to  say  if  there  is  a  single  best  receiver  synchronization 
time  t R  .  During  most  of  the  orbit  it  appears  that  the  time  tr  =  13.5  ms  would  be  the  best.  Flowever, 
this  value  of  tr  appears  to  give  a  value  of  |R,r(xR)|  that  is  severely  reduced  at  £  =  90°  and  £  =  270°. 
For  this  case  a  system  of  selection  of  three  values  of  tr  separated  by  1-ms  intervals  (12  5.  13  5,  and 
14.5  ms)  would  give  a  more  stable  average  performance  in  the  evaluation  of  R /(r*).  The  discon¬ 
tinuity  in  Fig.  93d  at  tr  =  13.8  ms  signified  by  the  heavy  vertical  line,  is  due  to  a  360°  ambiguity  that 
has  no  real  significance  except  at  the  nulls  (£  =  120°  and  £  =  240°),  where  iR/tr#)!  —  0.  Com¬ 
parison  of  Figs.  93(c)  and  93(d)  illustrate  this  feature  clearly. 

Figure  94  shows  the  correlation  performance  of  an  antenna  with  y  —  80c.  Near  th  =  14.0  ms 
and  £  =  180°  there  is  an  abrupt  drop  of  signal  and  jump  in  phase.  Otherwise  the  t R / ( r w )  ■  function 
seems  fairly  well  behaved. 

Figure  95  shows  the  behavior  of  R  r(rR)  lor  an  antenna  with  y  =  60°.  The  amplitude  and  phase 
variation  of  the  correlation  are  much  milder  than  for  the  cases  with  smaller  values  of  y,  since  the  verti¬ 
cal  electric-field  components  exceed  the  horizontal  components  by  at  least  one  order  of  magnitude  at 
this  droop  angle.  This  trend  is  continued  as  y  is  decreased  to  45°  and  then  to  0°  as  in  Figs.  96  and  97. 

Figure  98  shows  the  value  of  R  / 1  Kl'(r  K  )  It  illustrates  the  case  of  pure  first-order-mode  propaga¬ 
tion  from  a  vertical  electric-dipole  source. 

The  values  of  the  propagation  quantities  in  Table  7  are  similar  to  those  of  Table  6  except  for  an 
adjustment  to  produce  maximum  mode  cancellation  and  interference  effects.  We  have  made  the  first- 
order  and  second-order  vertically  generated  modes  have  opposite  phase  and  equal  magnitude  We  have 
modified  the  phase  of  the  horizontally  generated  modes  so  they  are  in  the  same  phase  plane  with  the 
vertically  generated  modes  to  produce  maximum  field-strength  variation  Hence  both  severe  effects 
appear  in  Figs.  99  through  104. 

Figure  99  gives  the  behavior  or  R  ’(tk)  for  the  propagation  parameters  of  Tabic  3  The  results 
arc  almost  identical  to  those  of  Fig.  92,  since  y  =  90°  and  the  only  changes  made  in  the  parameters  of  , 
the  horizontally  generated  modes  arc  to  the  phase  delay  times  r  and  rf1"'. 

Figure  100  shows  the  behavior  of  R/(t„)  for  the  case  of  y  =  85°.  The  influence  of  the  splitting 
of  the  vertically  generated  function  is  more  evident  on  the  phase  plot  than  on  the  amplitude  plot. 

Figure  101  shows  the  behavior  of  R/(tr)  for  the  case  of  y  =  80°  Here,  the  influence  of  the 
correlation-peak  splitting  of  the  vertically  generated  function  is  quite  evident  as  the  notch  in  the  corre¬ 
lation  peak  makes  its  way  from  side  to  side  through  the  values  of  rR  and  £.  This  is  illustrated  further 
in  Figs.  102  and  103,  which  give  the  graphs  for  y  “  60°  and  45°. 

Finally.  1  jg.  104  shows  the  split  correlation  peak  for  a  purely  vertical  antenna  being  carried  around 
an  orbit  Here  we  find  a  nearly  constant  variation  of  |R/(r,)|  versus  £,  and  the  phase  ol  the  correlation 
function  i i>i  (tk)  varies  only  moderately  with  £.  As  a  function  of  t# .  these  variables  show  the  expected 
splitting. 

DISTLRBANCF.S 

As  mentioned  previously,  long  wavelength  waves  require  an  ionospheric  reflection  medium  to 
provide  lor  their  world  ranging  capabilities.  This  reflecting  ionospheric  medium  is  caused  by  the 
natural  interaction  of  the  sun  and  the  stars  with  the  earth's  atmosphere,  mediated  by  the  interplanetary 
medium  and  by  the  geomagnetic  field  of  the  earth.  In  total  this  constitutes  an  extremely  complex  sys¬ 
tem  that  is  poorly  understood  and  monitored.  An  additional  complication  arises  in  the  consideration  of 
the  ground  —  r he  lower  boundary  of  the  waveguide.  The  effective  conductivity  of  the  ground  surface 
can  vary  over  great  ranges  The  conductivity  is  controlled  by  the  composition  and  layering  of  the  earth. 


104 


«! 


\N 


a 

sr% 


SS 


t: 

Nj 

s’ 


Ls 

,S 


I 

j 


I 

I 


NRL  REPORT  9028 


by  the  water  content  of  the  soil  and  rocks,  by  their  roughness,  and  by  their  temperatures  One  may 
consider  the  variation  of  ground  conductivity  with  location  as  a  disturbance  from  the  placid  intellectual 
picture  of  a  homogeneous  propagation  path.  Also,  the  seasonal  wetting  and  warming  of  vast  areas  of 
earth  should  be  of  consequence  to  radio  propagation  forecasts.  Further,  the  lower  atmosphere  contains 
elements  that  slightly  retard  the  propagation  of  radio  waves,  form  ducting  layers  and  account  for  the 
famous  "four  thirds  earth  radius”  law  that  becomes  applicable  even  at  low  frequencies.  Therefore,  the 
weather  is  an  additional  complicating  factor. 

Many  things  can  disturb  radio  waves,  but  to  what  degree  and  how  deleterious  the  various  distur¬ 
bances  are  is  not  completely  known.  There  are  some  answers,  however,  and  a  great  deal  of  investiga¬ 
tion  is  still  tequired. 

Solar  Flares— Sudden  Ionospheric  Disturbances 

Long-wavelength  waves  have  always  been  thought  of  as  more  reliable  than  high-frequency  waves, 
since  the  solar  flares  and  magnetic  storms,  which  tend  to  completely  obliterate  long-haul  high-frequency 
communications  circuits,  were  not  able  to  shut  down  communications  on  long  waves.  Figuratively,  the 
very  low  frequency  waves  effectively  “ducked  under"  the  enhanced  ionization  layer  of  the  lower  iono¬ 
spheres  that  cause  the  severe  high-frequency  signal  degradation  The  solar  flares  are  noticeable  on 
long-wavelength  circuits,  but  they  are  normally  seen  as  a  change  of  signal  level  and  not  as  a  complete 
loss  of  signal.  Sometimes  the  signal  will  increase ,  sometimes  it  will  decrease.  The  precise  behavior  is 
currently  unpredictable  without  making  exacting  predictions  of  the  locations  of  interference  nulls  under 
normal  and  disturbed  conditions  On  the  whole  I  believe  that  the  VLF  signal  is  more  likely  to  be 
enhanced  than  depressed  during  a  solar  flare  for  an  all-daytime  path  at  temperate  latitudes  [59], 

Polar  Cap  Absorption  Events 

Polar  cap  absorption  (PCA)  events  are  natural  disturbances  that  affect  VLF  propagation  for  ,'aths 
crossing  ihc  polar  region^.  PCAs  are  caused  by  a  llux  of  charged  particles  (protons)  emitted  from  the 
sun  during  a  solar  flare  The  particles  follow*the  lines  of  the  magnetic  field  and  effectively  lower  the 
height  of  the  ionosphere  in  the  polar  regions  because  of  the  enhanced  ionization.  The  onset  of  a  PCA 
effect  on  VLF  is  more  gradual  than  that  of  a  sudden  ionospheric  disturbance,  however  the  effects  are 
more  long-lasting  PCA  events  lasting  5  days  are  not  uncommon,  and  they  occur  fairly  frequently. 
There  were  1 1  PCAs  between  March  24,  1966  and  June  6,  1967. 

Papers  by  Wesierlund  el  al  (60|  and  by  Dalland  el  al.  [61 1  contain  detailed  discussions  of  PCA 
events  and  their  effects  or.  VLF  signals.  A  brief  summary  of  the  material  of  these  references  that 
seems  relevant  to  a  communications  channel  model  follows 

During  a  PCA  the  phase  of  a  VLF  signal  normally  advances  to  some  maximum  amount  and  later 
returns  to  its  pre-PC  A  value  Phase  advances  of  as  much  as  50  fi  s  have  been  observed.  The  behavior 
of  the  amplitude  of  a  VLF  signal  during  a  PCA  is  not  completely  predictable  For  paths  over  water  :n 
the  polar  regions  the  amplitude  variations  arc  generally  less  severe  and  unpredictable  in  sign.  However, 
for  paths  that  involve  a  good  deal  of  extremely  low  conductivity  soil,  Greenland  or  Antarctica,  the 
amplitude  consistently  shows  strong  signal  losses  (10  to  20  dB)  when  a  PCA  occurs  Figure  105  shows 
the  effects  ol  three  PCA  events  on  the  VLF  signal  over  three  Arctic  paths  Note  that  the  paths  NPG- 
Siockholm  and  WWVL-Stockholm  show  severe  attenuation  during  the  PCAs  Both  of  these  paths  pass 
over  the  low  conductivity  ground  of  Greenland  The  path  from  NPM-Stockholm  does  not  pass  over 
Greenland,  and  relatively  small  changes  in  amplitude  are  observed. 

The  behavior  of  the  V  LI  channel  during  PCA  event  is  veiy  significant  for  communications  over 
some  high-latitude  paths.  A  very  low  data  rale  modulation  scheme  could  compensate  for  the  very- 
severe  long  duration  attenuation  that  occurs  At  the  same  time  the  phase  variation  is  I'aiily  slow  over  a 
message  lime  and  should  not  pose  a  severe  obstacle  to  communication. 
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Fig  105  —  VI, F  phase  and  amplitude  data  for  the  transmissions  monitored 
during  PCAs  occurring  on  May  23,  1967,  May  28,  1967  3nd  June  6.  1967 


Transition  Fading 

Transition  lading  refers  to  signal  fades  that  occur  when  the  sunrise  or  sunset  line  (the  terminator) 
crosses  the  propagation  path  between  the  transmitter  and  the  receiver  This  phenomenon  has  been 
studied  by  Crombie  [62-64]  and  others  [65. 66).  We  summarize  here  the  material  from  these  refer¬ 
ences. 

In  Crombie’s  model,  the  difference  in  ionospheric  heights  on  either  side  of  the  terminator  causes 
modal  conversion  in  the  waveguide.  The  resultant  interference  on  the  side  of  the  terminator  away  from 
the  transmitter  gives  rise  to  the  transition  fading.  The  principal  effect  of  the  fading  is  to  introduce  mul¬ 
timode  complications  similar  to  those  already  observed  at  nighttime  During  these  times  the  amount  of 
the  multimode  interference  is  determined  by  the  relative  amplitudes  of  the  modes  and  the  group  and 
phase  delays  involved.  The  effect  should  be  comparable  to  that  observed  during  normal  nighttime  con¬ 
ditions  wh  n  mullimode  effects  are  also  important. 
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The  transition  fading  is  a  dynamic  phenomenon  because  the  terminator  moves  along  the  path 
length.  Therefore,  the  rales  of  change  of  phase  and  amplitude  are  important.  Crombie  [62]  gives  an 
observed  maximum  dd>ldt  of  —70  ±27min.  However,  such  phase  rates  could  become  very  large  for 
very  short  times,  and  a  step  function  jump  of  phase  is  possible  if  a  very  sharp  and  severe  modal 
interference  null  occurs  167],  The  rate  of  change  of  amplitude  with  time  is  approximately  0.3  dB/min 
during  a  fade.  Faster  rates  of  change  are,  of  course,  possible  for  very  deep  fades.  Typical  fade  depth  is 
approximately  10  dB;  however,  greater  fades  have  been  observed,  and,  of  course,  even  perfect  cancella¬ 
tion  of  modes  is  theoretically  possible. 

Schneider  and  Carlston  [66]  performed  a  statistical  analysis  of  transition  fading  depths  and  dura¬ 
tions.  They  have  applied  the  statistics  of  extreme  values  [68]  to  data  of  Rhoads  and  Garner  (69)  for  the 
path  from  Haiku  to  Washington,  DC,  that  occurred  in  May  and  June  1965.  They  found  a  good  fit  to  a 
probability  distribution  for  depth  of  fade  and  outage  times  given  by  type  I  and  type  HI  distributions  of 
the  statistical  theory  of  extreme  values 

The  cumulative  probability  distributions  for  the  depths  of  fade  and  the  times  of  signal  strength,  10 
dB  below  the  normal  daytime  level,  are  shown  in  Figs.  106  and  107.  Experimental  examples  of  the 
diurnal  variation  of  received  signal  from  Haiku  to  Washington,  DC  [69],  are  shown  in  Figs.  108  and 
109. 
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FRACTION  OF  FADE  MINIMA  GREATER  THAN  ORDINATE 

Kig  I  OP  —  Distribution  of  transition  fades  for  10  consecutive  days, 
distribution  of  the  third  sunrise  fade  for  IV  consecutive  days 
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Fig  109  —  NPM  (26  1  kc/s)  data  recorded  on  May  21,  22.  23,  24, 
near  Washington,  DC 

Electron  Precipitation-Induced  Perturbations 

Anyone  who  has  recorded  longwave  field  strengths  will  note  that  the  nighttime  traces  are  more 
irregular  and  "shakey"  than  the  daytime  ones  — even  though  the  Helds  are  usually  stronger  and  less 
affected  by  ambient  noise  [70!. 

The  picture  that  is  emerging  to  explain  this  feature  [71]  of  propagation  behavior  is  curious.  It 
appears  that  the  earth's  magnetosphere  is  filled  with  energetic  particles  in  the  famous  Van  Allen  belts. 
These  belts  are  continually  filling  up  with  particles.  The  equilibrium  reached  depends  on  the  plasma 
density  in  the  magnetosphere  and  the  intensity  of  lightning  generated  whistler  waves  that  enter  the 
magnctospheie  from  the  earth  ionosphere  waveguide  [72,73).  The  propagation  of  whistler  waves 
through  the  magnetosphere  can  shake  or  scatter  the  trapped  electrons  into  their  "loss  cone,"  and  they 
piecipilate  from  the  Van  Allen  Belts.  When  they  fall  from  the  magnetosphere,  they  strike  the  lower 
atmosphere  and  produce  electrons  in  the  nighttime  D-iayer.  This  sudden  electron  production  causes  a 
quick  change  in  the  ionospheric  reflection  height  and  thereby  perturbs  the  waves  propagating  in  the 
earth-ionosphere  waveguide  causing  the  abrupt  events  noted  in  nighttime  recordings. 

The  reader  is  invited  to  learn  more  about  this  "strange  but  (probably)  true"  mechanism  by  study¬ 
ing  Refs.  74  through  79. 

The  auroral  zone  is  a  region  where  the  electron  precipitation  is  greatly  intensified  over  midluti- 
tudes.  It  should  be  expected  that  these  electrons  will  produce  disturbance  also  (70,80-83). 
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ELF  Nighttime  and  PCA  Disturbances 

Tiie  ELF  longwave  communications  system  has  an  interesting  disturbance  mode  that  has  not  been 
exhibited  by  the  VLF  waves.  Because  of  the  huge  wavelength  of  ELF  (5-3000  Hz)  waves,  their  reflec¬ 
tion  coefficient  matrix  (Eq.  (D)  at  nighttime  depends  strongly  on  the  electron  density  distribution  up 
to  the  F-layer!  This  means  that  the  ELF  nighttime  wave  will  be  subject  to  the  vagaries  that  HF  exhibits 
[84,85]  —  except  the  ELF  will  be  averaged  over  a  larger  Fresnel  zone.  This  possibility  was  pointed  out 
early  by  J.  Galejs  [86-88!.  Now  in  recent  years,  it  is  being  treated  in  greater  detail  and  with  greater 
respect  [89-91]. 

Reference  92  treats  the  interesting  possibility  that  the  change  in  ELF  propagation  velocity  under  a 
strong  PCA  could  bend  the  ELF  waves  so  as  to  cause  weakening  and  reenforcements  of  the  waves  in 
unpredicted  ways. 

C-Layer  Disturbance 

Under  some  circumstances  during  daytime,  cosmic  rays  can  produce  a  significant  ionization  layer 
beneath  the  D-layer.  This  is  known  as  a  C-layer  [93.94],  The  formation  and  disappearance  of  this 
layer  is  somewhat  mysterious  (but  so  is  the  D-layer).  Obviously,  with  the  right  intensity  this  layer  will 
cause  its  own  assembly  of  effects.  The  most  notable  effect  is  a  predicted  extra  attenuation  for  waves 
that  go  twice  through  this  layer  before  reflecting  at  the  D-layer  Cosmic  rays  are  known  to  vary  in 
intensity  with  the  11 -year  solar  cycle  and  with  solar  eruptions  [95]. 

Ground  Conductivity  and  Weather  Effects  in  Low-Frequency  Radio  Propagation 

In  the  low-frequency  propagation  research  community  awareness  is  growing  of  the  need  for  a 
better  propagation  model  that  accounts  for  ground  conductivity,  terrain,  and  atmospheric  conditions 
[96-102].  For  example.  Fig.  1 10  from  Ref.  99  shows  the  strong  relationship  between  the  time  of  arrival 
of  a  Loran-C  groundwavc  signal  transmitted  from  Carolina  Beach,  NC,  to  Fort  Wayne,  IN.  and  the 
atmospheric  temperature  recorded  at  Roanoke,  VA,  very  close  to  this  propagation  path.  Reference  96 
puts  forth  theoretical  ideas  to  explain  the  weather  related  propagation  delay  in  terms  of  changes  in  the 
gradient  of  atmospheric  index  of  refraction.  The  disturbance  of  radio  wavefields  caused  by  irregular 
terrain  and  ground  conductivity  variations  is  often  observed.  For  example.  Fig.  Ill  from  Ref.  103 
shows  data  illustrating  the  "recovery  effect"  observed  in  passage  over  ground  segments  that  have  dis¬ 
similar  conductivities.  Figure  111(a)  shows  the  flight  path  of  a  receiving  aircraft  over  Great  Britain. 
Figure  111(b)  shows  the  dramatic  "recovery  effect"  as  the  groundwavc  passes  over  a  river  and  over  a 
land/sea  boundary  The  experimental  path  was  chosen  because  of  the  flatness  and  homogemety  of  the 
terrain,  otherwise  even  greater  irregularities  would  be  evident  in  the  portion  of  the  flight  over  land 

Figure  112  1104]  illustrates  the  variation  of  atmospheric  index  of  refraction  vs  height  which  can 
produce  ducting  at  radar  frequencies.  Such  atmospheric  changes  should  produce  measurable  effects  at 
the  edge  of  the  effective  coverage  range  of  a  groundwave  communications  or  navigation  system. 
Indeed  the  100  kHz.,  Loran-C  community  has  noticed  this  effect  (105). 

NOISE 

The  prediction  of  atmospheric  noise  is  a  very  important  portion  of  the  overall  prediction  of  long 
wavelength  system  performance  The  basic  source  of  most  long  wavelength  system  noise  is  lightning  in 
thunderstorms  During  the  1960s,  the  National  Bureau  of  Standards  conducted  a  long  scries  of  atmos¬ 
pheric  noise  measurements  at  stations  around  the  world  that  led  to  the  CC'IR  322  publication  1106,107). 

Ma  xwell  ard  others  of  the  Westinghouse  Georesearch  Corporation  utilized  NBS  data  along  wit’,, 
propagation  model,  and  the  statistical  records  of  thunderstorm  occurrence  to  make  a  more  realistic 
model  of  noise  for  toe  frequency  range  from  10  to  30  kHz  1108)  This  model  was  further  improved 
and  extended  to  predict  atmospheric  noise  ■!  any  altitude  and  polarization  within  the  earth-ionosphere 
waveguide  and  also  up  to  nearly  60  kHz.  (109-1  12). 
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The  NBS  noise  data  (106)  was  taken  through  the  LF  and  medium  frequencies  up  to  the  HF  range. 

At  ELF  frequencies,  similar  studies  have  been  made  (114-119)  on  a  sporadic  basis.  No  official  group 
has  released  worldwide  maps  of  ELF  noise  and/or  powerline  interference.  It  would  seem  appropriate  to 
extend  the  model  of  Ref.  112  down  to  ELF  and  up  through  MF  to  HF  frequencies. 

Some  proportion  of  the  noise  is  generated  by  currents  flowing  in  the  ionosphere  and  magneto¬ 
spheres  (120).  These  noise  sources  are  believed  to  be  more  important  at  polar  receiver  sites  throughout 
the  longwave  band,  but  a  clear  understanding  of  precisely  when,  where,  and  how  much  noise  comes 
from  extraterrestrial  sources  is  not  known  (113). 

Noise  Mitigation 

Because  the  highest  noise  pulses  generated  by  lightning  are  strongly  nonthermal  in  amplitude  pro¬ 
bability  distribution  various  methods  have  been  used  to  reject  the  energy  contained  in  these  pulses  so 
as  to  reduce  the  “effective”  root-mean-square  (RMS)  noise  in  a  receiver  much  below  that  which  would 
be  seen  by  a  “true”  RMS  noise  meter.  It  is  believed  that  Shannon’s  Theorem  on  Communication 
Channel  Capacity  is  valid  for  gaussian  noise  statistics  but  that  for  non-Gaussian  noise  receiver  can  do 
better  than  look  at  the  RMS  noise.  Noise  theorists  have  defined  a  parameter  called  Vd  that  is  the  ratio 
of  the  RMS  noise  envelope  to  the  average  noise  envelope.  It  is  believed  that  an  effective  noise  clipper 
circuit  can  reduce  the  noise  level  down  to  about  the  average  levels.  (This  belief  is  always  being  sorely 
tested  by  the  untractable  performance  of  real  systems  in  the  real  noise  world.)  Caution  is  advised  for 
anyone  who  thinks  that  a  few  sim^.c  measurements  will  suffice  to  characterize  a  receiver's  performance 
in  noise  for  all  time.  Noise  is  noisy!  It  is  variable  in  time  and  place,  and  a  good  amount  of  experience 
with  any  nonlinear  receiver  must  be  obtained  before  its  performance  can  be  relied  on.  It  is  even  ques¬ 
tionable  how  io  characterize  a  receiver's  performance  in  noise.  Noise  has  an  amplitude  probability  dis¬ 
tribution  and  a  time  probability  distribution;  but  these  distributions  may  not  be  enough  to  characterize 
noise  to  permit  a  receiver’s  performance  to  be  specified. 

SYSTEM  PERFORMANCE  CALCULATION 

At  ELF,  Ref.  121  gives  a  good  example  of  a  longrange  communication  system  performance  calcu¬ 
lation  At  VLF  and  LF,  the  system  performance  calculation  is  facilitated  by  computer  programs 
(38,46,47)  that  can  be  used  to  produce  coverage  maps  such  as  those  illustrated  in  Figs.  56  and  57  and 
in  Figs.  113  through  116.  One  merely  decides  what  signal  to  noise  ratio  (S/N)  is  needed  for  reliable 
communication  for  a  given  keying  rale  receiver  and  antenna.  One  titer  selects  the  appropriate  contours 
that  describe  this  coverage. 
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Fig.  113  —  Signal  level 
contours  in  dB  above  I 
liV/m  for  the  NAA 
transmitter  in  Cutler,  ME 
(Lai  44  60°,  Long 
-67  20°)  for  I  MW 
radiated  power  during  July 
at  24  0  kHz  jnd  90%  lime 
a\  ailability 


Fig.  114  —  Signal  level 
contours  in  dB  above  I 
jxV/m  for  the  NAA 
transmitter  in  Cutler,  ML 
(Lat  44,t>0°,  Long 
-67.20°)  for  1  MU 
radiated  power  during  July 
at  24.0  kHz  and  99 Vo  time 
availability 
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Appendix 

THE  IMPORTANCE  OF  THE  ELF/VLF/LF  BAND 


The  Ground  Wave  Emergency  Network  (GWEN)  will  provide  a  nuclear  hardened  terrestrial  com¬ 
munications  system  to  carry  warning  data  and  emergency  action  messages  and  force  direction  orders  to 
U.S.  forces.  When  the  first  phase  is  operational,  in  FY  1987,  GWEN  will  be  a  unique,  low-cost,  har¬ 
dened  communications  system  serving  both  national  and  nuclear  CINC  requirements.  The  second 
phase  of  the  program  will  have  additional  transmitters  and  receivers,  and  will  provide  a  hardened  link  to 
the  intercontinental  ballistics  missile  force  and  other  strategic  elements.  We  expect  to  procure  about  50 
more  nodes  for  the  network  in  FY  1987  and  to  complete  the  system  by  the  end  of  the  decade. 

To  improve  communications  connectivity  with  strategic  bombers  in  flight,  we  plan  to  equip  them 
with  Miniature  Receive  Terminals  (MRTs)  for  receiving  low-frequency  (LF)  and  very  low  frequency 
(VLF)  communications.  Although  slower  than  transmissions  in  higher  frequency  bands,  LF/VLF  com¬ 
munications  can  be  transmitted  over  much  greater  distances  than  can  line-of-sight  communications,  and 
they  are  much  less  susceptible  to  nuclear  effects  and  jamming  than  the  existing  satellite  and  high- 
frequency  (HF)  communications  are.  Flight  tests  of  the  terminals  are  planned  to  be  completed  in  FY 
1987,  production  should  be  completed  in  early  FY  1988. 

To  communicate  with  ballistic  missile  submarines,  the  Navy  maintains  special  radio  relay  aircraft 
continuously  airborne  over  the  Atlantic  and  Pacific  oceans.  These  planes,  called  Take  Charge  and 
Move  Out  (T  AC  AMO),  transmit  relayed  messages  to  the  submarines  by  a  variety  of  communications 
media.  In  FY  1989,  we  will  begin  replacing  the  EC-130s  (modified  C-130  transports)  now  flying  the 
TACAMO  mission  with  the  faster,  more  survivable,  more  endurable  and  longer  range  E-6A  (a  deriva¬ 
tive  of  the  Boeing  707  airborne  warning  and  control  system  I  AW  ACS]  airframe).  At  the  same  time, 
the  SSBN  force  will  be  outfitted  with  improved  VLF  receivers.  The  E-6  A’s  ability  to  operate  over  vast 
expanses  of  the  ocean  wili  permit  SSBNs  and  nuclear  cruise  missile  equipped  submarines  to  expand 
their  operating  areas  and  still  receive  messages  from  the  President.  Building  toward  a  fleet  of  15  E- 
6As,  we  are  requesting  funds  to  procure  the  third,  fourth,  and  fifth  production  aircraft  in  FY  1987. 
Although  the  planes  will  be  outfitted  initially  with  communications  equipment  transferred  from  EC- 
130s,  they  eventually  will  carry  a  vastly  improved  VLF  communications  system. 

Submerged  submarines  now  must  deploy  an  antenna  at  or  close  to  the  ocean’s  surfa'  ?  to  receive 
messages,  thus  potentially  increasing  their  susceptibility  to  detection.  The  use  of  extrerr  Jy  low-  fre¬ 
quency  (F.LF)  communications,  which  can  penetrate  to  great  depths,  will  alleviate  this  constraint.  The 
FLF  Communications  System  consists  of  two  transmitlei  sites  (in  Wisconsin  and  northern  Michigan), 
operating  in  electrical  synchronization  with  receivers  aboard  submarines.  The  ELF  system’s  high  relia¬ 
bility  and  continuous  transmissions  significantly  will  upgrade  peacetime  communications  to  deployed 
submarines  and  will  support  the  transition  to  wartime  operations.  Receivers  are  being  installed  in  some 
submarines  this  year  so  that  the  full  system  can  be  tested  in  FY  1987.  We  are  requesting  funds  for  FY 
1987  to  procure  receivers  for  other  submarines.  The  system  should  be  fully  operational  by  the  end  of 
FY  1988 
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